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INTRODUCTION 


This  report  is  a  continuation  of  the  author's  "Theories  of 
Nutation  and  Polar  Motion  I",  Report  No.  309,  Department  of 
Geodetic  Science,  The  Ohio  Stabe  University,  1980.  It  considers 
theories  based  on  earth  models  consisting  of  an  elastic  mantle 
and  a  1 iquid  core . 

The  first  theoretical  treatment  of  such  a  model  has  been 
given  already  by  Poincarfc  (1910),  but  only  to  the  extent  to  show 
that  a  nearly  diurnal  free  wobble  (NDFW)  exists.  His  treatment 
is  sketchy  and  incomplete;  his  approach  has  not  been  followed  by 
later  authors.  No  quantitative  results  are  given. 

Only  in  1949,  Harold  Jeffreys  took  up  the  problem  again  and 
used  a  variational  method,  although  a  somewhat  simplified  one, 
to  get  preliminary  numerical  values.  A  perfected  version  of  the 
variational  approach  (Jeffreys  ar<d  Vicente,  1957a, b)  gave  the  first 
numerical  results  that  are  so  detailed  and  accurate  as  to  come 
close  to  modern  standards,  although  the  underlying  earth  models 
were  rather  simplified. 

The  variational  approach  of  Jeffreys  and  Vicente  has  not  been 
pursued  since.  Molodensky  (1961)  writes:  "Jeffreys  and  Vicente 
applied  a  variational  method  which  made  it  difficult  to  have  a 
clear  notion  about  the  degree  of  approximation  to  the  exact  so¬ 
lution  of  the  problem.  I  attempted  to  reproduce  this  theory  and 
to  make  it  more  exact,  but  owing  to  numerous  obscurities  I  had 
to  give  up  this  method." 

Thus  Molodensky  (1961)  developed  a  method  of  his  own,  based 


v 


on  a  spherical  approximation  of  the  equations  of  elasticity  for 
the  mantle  (and,  possibly,  for  a  solid  inner  core)  and  on  a 
hydrodynami cal  theory  for  the  liquid  core.  He  got  numerical  re¬ 
sults  which  are  excellent  even  by  present  standards;  furthermore, 
this  approach  --  generalized,  modified,  simplified  --  has  been 
the  starting  point  of  practically  all  subsequent  work. 

Unfortunately,  the  classical  papers  by  Jeff reys-Vi cente 
and  Molodensky  are  very  difficult  to  read,  because  of  "certain 
expository  difficulties  in  both",  as  Smith  (1977)  remarks.  On 
the  other  hand,  they  are  important  enough  to  deserve  careful 
study. 

The  present  report  attempts  an  elementary  introduction  to 
the  work  of  Jeffreys,  Molodensky  and  subsequent  developments.  We 
shall  try  to  interrelate  various  approaches  by  using  a  basic  set 
cf  parameters  (essentially  Mol odensky ' s )  and  relating  to  it 
other  variables  such  as  employed  by  Jeffreys.  The  various  co¬ 
ordinate  systems  used  by  different  authors  will  be  seen  to  redu¬ 
ce  to  two  basic  systems,  the  nutation  frame  (Jeffreys,  Smith, 
Wahr)  and  the  body  frame  (Molodensky,  Shen-Mans i nha ) ,  which  are 
easily  related. 

The  most  difficult  part  of  (Molodensky,  1961)  is  the  hydro- 
dynamic  theory  of  the  core.  In  order  to  clarify  the  essential 
features  and  to  supply  motivations  which  are  missing  in  Molo- 
densky's  presentation,  we  present  a  pedestrian's  version  using 
a  simplified  model  of  a  homogeneous  core.  This  model,  although 
inadequate  from  a  practical  point  of  view,  is  best  suited  for 
a  first  approach;  it  permits  an  instructive  comparison  of  the 
methods  of  Jeffreys  and  Molodensky  with  the  classical  theory  of 
Poincare  (for  a  rigid  mantle).  We  shall  also  show  the  relation 
of  Mol odensky ' s  theory  to  modern  treatments  of  the  core  in  terms 
of  toroidal  oscillations  (Smith,  Shen  and  Mansinha,  Wahr). 
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PART  A 


HYDRODYNAMICS  OF  A  HOMOGENEOUS  CORE 


1 .  Hydrodynamic  Equations 


In  an  inertial  frame  XYZ,  the  motion  of  an  ideal  (non- 
viscous) fl uid  is  described  by  the  well-known  equations 

X  «  gradV  -  pmlgradp  ;  (1-1) 

cf.  (Lamb,  1932,  pp.  3-4;  lanczos,  1970,  p.  36:.)-  Here 


X  -  d2X/dt2  (1-2) 

denotes  the  second  derivative  of  the  position  vector  X  *  (X,Y,Z) 
with  respect  to  time;  it  represents  the  acceleration  of  a  moving 
particle,  d/dt  denoting  a  Lagrangian  differentiation  following 
t!.e  motion  of  the  particle  (as  opposed  to  an  Eulerian  differentia' 
tion  referred  to  a  poi  ,it  of  fixed  coordinates).  The  vector  gradV 
is  the  gradient  of  a  potential  V  , 
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grad  V  «  (av/ax,  aV/aY,  av/aZ)  ,  (1-3) 

p  is  the  density  and  p  the  pressure  of  the  fluid. 

In  a  moving  frame  xyz,  rotating  with  an  angular  velocity 
vector  u,  equations  (1-1)  are  to  be  replaced  by 

••  •  * 

x  +  +  2w* )c  +  u* (u* x.)  ■  grad  V  -  p“lgradp  ,  (1-4) 

where  on  the  left-hand  side  there  occur  additional  accelerations, 
which  may  be  considered  apparent  or  inertial  forces: 


wxx  .  Euler  force, 

« 

2w*x  .  Coriolis  force,  (1-5) 

w*(wxx)  .  centrifugal  force; 


cf.  (Lanczos,  1970,  p.  101).  The  cross  (*)  denotes  the  vector 
product  as  usual . 

Let  us  take  the  z-axis  coinciding  with  the  rotation  axis; 
then  u  has  the  form 

«  -  (0,0, a)  .  (1-6) 


We  furthermore  assume  constant  rotation  so  that  w  ■  0.  Then 
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(1-4)  takes  the  form 


The  vector  (n2x,  n2y,  0)  Is  the  gradient  vector  of  the  centri 
fugal  potential 

♦  *  ^n2(x2+y2  )  .  (1-8) 

We  furthermore  assume  a  homogeneous  fluid  with 

p  ■  const.  ,  (1-9) 

and  we  introduce  the  function 

*  V  +  ♦  -  p”*p  (1-10) 

and  the  displacements 
u  «  x  -  x0  » 

v  «  y  -  y0  »  u  »  £  -  xo  (1~H) 

W  ■  Z  -  Zq  , 


4 


of  a  fluid  particle  with  respect  to  some  fixed  initial  position 
x0  so  that 

•  *  II  M 

£  ■  JU  ,  X  ■  _U  .  (1*12) 


Then  the  system  (1-7)  reduces  to 


••  • 

U  -  2QV  ■  -3iJt/3X  , 

••  * 

v  +  2flu  ■  -3iji/3y  , 

•• 

W  *  -  3ij»/  32  , 

or  in  vector  notation, 

♦*  * 

u  +  2^x_u  «  -  grad  ♦  . 


(1-13) 


(1-14) 


This  form  has  been  used  by  Hough  (1Pj5,  p.  479)  in  his  pioneering 
investigation  of  an  earth  model  featuring  a  rigid  mantle  and  a 
liquid  core;  cf.  also  (Moritz,  1980b, p.  100). 

An  important  special  solution  of  the  system  (1-13)  is  ob¬ 
tained  by  putting 


<1*  ■  B(xz  cos  ot  +  yzsinat) 


(1-15) 


with  constants  8  and  o  .  Then 
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B 

u  -  -  z  cos  t 

<3(0+ 2a) 


v  - - z  s i r>  at  , 

o ( a+2P ) 


w  *  —  (x  cos  at  +  y  sin  at) 
a2 


(1-16) 


satisfy  (1-13),  as  we  immediately  verify  by  substitution. 

Nutation  frame  and  body  frame.  In  the  present  report,  we 
shall  use  two  basic  reference  frames  in  which  the  earth's  mantle 
is  approximately  at  rest.  The  first  is  the  body  frame  x!x2x3 
which,  for  a  rigid  mantle,  is  rigidly  connected  to  it  in  such  a 
way  that  the  x3  axii  s  the  symmetry  axis  (figure  axis)  of  the 
earth,  supposed  to  be  an  ellipsoid  of  revolution.  For  an  elastic 
mantle,  x3  represents  the  Liouville  axis  which  corresponds  to 
the  figure  axis  of  the  undeformed  earth;  it  is  the  same  as  the 
z-axis  in  (Moritz,  1980b,  p.  49).  This  frame  is  appropriate  for 
polar  motion,  which  is  motion  around  the  x3-axis. 

In  (Moritz,  1980b,  pp.  76)  we  have  introduced,  in  addition 
to  the  body  frame  xix2x3  ,  an  auxiliary  system  x°x2x°  which 
is  connected  to  the  inertial  system  in  a  prescribed  way:  the  x3- 
axis  has  a  fixed  direction  in  inertial  space,  and  the  system 
x?x°2x?  rotates  with  constant  angular  velocity  n  around  the  x3- 
axis.  This  system  is  appropriate  for  the  treatment  of  nutation. 


".'fa.--:,...  .V-. 


HffPBHW 


uriPgL'  2 
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which  is  motion  around  the  x^axis  (Moritz,  1980b,  p.  92).  There 
fore,  the  auxiliary  uniformly  rotating  system  x°x2X3  will  be 
called  nutation  frame. 

The  body  frame  and  the  nutation  frame  differ  only  by  a 
small  rotation: 


x  *  (I  +  £)x° 


where 


0 

■0  3 

e2 


93 

0 

"0 1 


-Q. 


(1-17) 


(1>18) 


is  a  small  skew-symmetric  matrix. 

In  view  of  (1-6)  it  is  clear  that  the  equation  of  motion 
in  the  nutation  frame  has  the  form  (1-14). 

Equation  of  motion  in  the  body  frame.  This  equation  has 
a  different  form.  To  derive  it,  we  put 

u  =  w0  +  ,  ( 1 - 1 9  a ) 

w o  *  (0,  0,  n)  ,  *  n(mi,  m2  ,  m3)  .  ( 1  - 1 9  b ) 


.v-rWtW'< 


— tisssBK, 1  •■ 
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We  shall  only  deal  with  polar  motion,  which  is  the  change  of 
di recti  on  of  the  vector  w  described  by  mi  and  m2,  but  not 
with  i rregul ari ties  in  the  angular  velocity,  which  is  the  change 
of  magni tude  of  the  vector  w  described  by  m3.  In  agreement 
with  (Moritz,  1980b, pp.  19  and  85),  we,  therefore,  put 


m3  =  0  .  (l-19c) 

The  use  of  the  new  system  xyz  means  that  we  have  to  con¬ 
sider  5u  and  add  its  effect  to  the  left  hand  side  of  (1-7)  which 
by  (1-6)  corresponds  to  w  *  u0  .  By  (1-5),  we  have  an  Euler 
force 


*  6uxX  , 


(1-20) 


a  change  of  Coriolis  force,  using  (1-12), 

•  • 

6  f^0  =  2  <$  w  *  X.  *  2  6(£X U_  ,  (1-21) 

and  a  change  of  centrifugal  force. 


*  -* x-)  “  «ox(«oxx)  •  (1-22) 

We  shall  consider  6  uj_  and  u_,  as  well  as  their  derivatives, 
as  quantities  of  first  order  and  shall  consistently  neglect  terms 


of  second  and  higher  order.  Then  (1-21),  being  a  product  of 
first  order  terms,  is  of  second  order  and  will  be  neglected. 
Thus , 


4*co 


0  . 


To  the  same  approximation. 


(1-23) 


6f^  *  6ojx(u)o*x)  +  wqx($w*x)  ,  (1-24) 

which  is  found  by  substi  tuti  ng  ( 1- 13a )  into  (1-22)  and  neglecting 
the  quadratic  term  6u*(6u*x)  . 

We  work  out  the  vector  products  in  formulas  (1-20)  and  (1-24) 

using 


*  (Qm^  ,  flm2  »  0) 


d-25) 


by  (l-19a,b).  The  result  is 


z 


5£c 


n2m2z 

n2(mix  +  m2y) 


d-27) 


3 

I 

-i 


(We  shall  indifferently  use  row  or  column  notation  for  vectors, 
whichever  is  more  convenient.) 

We  put 


IT^  *  eCOSot 

m2  »  csinat 


d-28) 


Since  and  m2  are  the  rectangular  coordinates  of  the  pole 

of  rotation  in  the  tangential  plane  at  the  unit  sphere  at  its 
intersection  with  the  z-axis,  (1-28)  corresponds  to  a  circular 
polar  motion  with  (angular)  frequency  o  and  radius  e  .  The 
situation  is  comparable  to  Eulerian  motion  illustrated  in  (Moritz 
1980b,  p.  11,  Fig.  2.1),  but  e  and  a  are  general  and  different 
from  a  and  og 

Then  (1-27)  becomes 


n2c  z  cosot 
n2e  z  si  not 

fl2e  ( x  cosot  +  y  si  n  o  t ) 


(1-29) 
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which  can  be  expressed  in  the  form 


=  -  grad  +  , 


(1-30) 


where 

i 


<t>  ■  -n2e(x2  cos  ot  ♦  yzsinot)  (1-31) 

is  an  incremental  centrifugal  potential.  The  Euler  force  (1-26) 
on  using  (1-28),  takes  the  form 


floe  Z  COS  ot 
fioeZ  sin  at 

-nae(xcosot  +  ysinot) 


(1-32) 


which  can  be  expressed  as 


i,  •  -  5  +  2  lUiJ  •  d-33) 

where  £3  »  (0,0,1)  is  the  unit  vector  of  the  z-axis.  It  is  im 
mediately  verified  that  (1-33)..  with  (1-31),  gives  (1-32). 

Thus  the  quantity 


n 


-e  +  6lc  c  1  +  9rad*  +  2  ir  II  -3  ( 1-34) 

must  be  added  to  the  left  hand  side  of  (1-14),  in  which  w  is 
understood  as  «o .  This  gives 

i 

j±  +  2u0*u  -  ( 1  +  2)  grad$  +  *  ft  ’H'  ®.3  *  -grad$o  »  ( 1-35) 

where  $  of  eq .  (1-10)  is  now  denoted  by  $o  : 

-$0  *  V  +  *  -  p_1p  ;  (1-36) 

*  is  again  defined  by  (1-8).  Let  us  finally  define  a  new  function 

*  by 

- V  *  -<l»o  +  (1+jjH  ■  V  +  *  -  p_1p  +  (1+2)$  .  '1-37) 

Then  (1-35)  takes  the  final  form 

u  +  2u0*u  »  -grad$  -  2  2  |2  e3  ,  (1-38) 

which  will  be  used  in  the  following  section. 


12 


2 .  Method  of  Molodens^y 


Molodersky's  method  is  not  restricted  to  fluids  of  constant 
density;  for  a  homogeneous  core,  however,  it  becomes  essentially 
simpler.  Thus  we  shall  here,  and  throughout  the  present  report, 
presuppose  a  core  consisting  of  an  ideal  homogeneous  fluid.  This 
will  provide  an  easy  introduction  to  an  otherwise  very  difficult 
theory  and  thus  facilitate  the  reading  of  presentations  such  as 
(Molodensky,  1961;  Oobert,  1964;  Melchior,  1978,  sec.  6.2). 

The  system  of  differential  equations  of  fluid  motion  (1-38) 
may  be  written 


u  -  2nv  - 
v  +  2nu  « 

w  » 


4*  • 

«y 

3  2  °  3  $ 

32  "  n  32 


where 


-i>  *  V  +  *  -  p_lp  +  ( 1  +  £)♦  , 

♦  ■  +  y2)  , 

«  -ft2  e ( xz  cosot  y  z  sinot)  , 


(2-1) 


(2-2) 

(2-3) 

(2-4) 


■BBBB 
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V  is  the  potential  of  gravitational  attraction  acting  on  the  fluid 
module  under  consideration,  and  p  is  the  pressure;  the  density  p 
is  constant. 

The  coordinate  system  xyz  is  such  that  the  mantle,  on  the 
average,  is  at  rest  with  respect  to  it,  and  the  z-axi$  has  a  con¬ 
stant  direction  with  regard  to  inertial  space.  The  unit  vector  of 
the  instantaneous  rotation  axis  has  the  x  and  y  coordinates 


mi  «  ccosot  ,  m2  ■  esinot  .  (2-5) 

If  the  rotation  axis  coincides  with  the  2-axis,  for  c  *  0  ,  then 
eqs.  (2-1)  reduce  to  (1-13). 

Let  us  now  consider  harmonic  oscillations  of  the  form 


u^(x,y,z,t)  «  £1  (x,y ,z)cosot  +  J±2(*»y ,z)sinot  ,  (2-6) 


depending  on  time  with  (angular)  frequency  0  .  Then 

u  -  -o2u  ,  ( 2-7 ) 

that  is, 

u  *  -0 2 u  ,  v  *  -o  2  v  ,  w  *  -a2w  ,  (2-8) 

as  we  directly  see  by  substitution.  Thus  (2-1)  becomes 
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<j2u  +  2nv  *  <|*  , 

CV 

o2v  -  2flu  ■  n,  ,  (2-9) 

«2»  ■  *,  +  2  §  ♦,  • 

D1 f ferenti itlon  with  respect  to  time,  using  (2-8),  gives 

•  • 
o2u  -  2flo2v  *  , 

.x  (2-10) 

o2v  +  2flo2u  *  \|i  , 

y 

where 

**  •  lx  *  *x  *  *  etc-  <2“11) 

•  • 

We  solve  (2-xn)  for  u  and  v  and  substitute  in  (2-9).  This 
gives  us 

(o2  -  4o2)u  *  *x  ’  2  , 

Vo2  -  402)v  «  *y  +  2  ^7  ,  (2-12) 

°2"  -  ♦«  *  2  ♦,  • 
which  solves  (2-1)  for  u,  v,  w  . 


15 


The  volume  dilatation  is  defined  by 


3u  .  3V  .  9W 
7x  +  TyTz 


(2*13) 


it  characterizes  the  change  of  a  volume  element  of  the  fluid  during 
motion.  A  homogeneous  fluid  for  which 


p  •  const. 


(2-14) 


in  space  and  time,  is  necessarily  Incompressible.  This  means 

e  •  0  .  (2-15) 

We  take  u,  v,  w  from  (2-12),  substitute  into  (2-13)  and  put  the 
result  equal  to  zero: 


-  2  *j»  +  ip  +2  )  + 

a2  -  4Q2  **  o2  *y  yy  a 2  xy; 


(*, 


_Lu  +  2  £  $  )  ■  0  . 

2  n  ZZJ 


Since  ■  0  by  (2-4),  we  are  left  with  Poincarfe's  equation 


2 

tp  +  *  +  (1-4  ■  0 

xx  yy  1  ,]  p ;  zz 


(2-16) 
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(Poincare,  1885;  Hough,  1895;  Jeffreys,  1949;  Molodensky,  1961). 

This  equation  is  similar  to  Laplace's  equation,  to  which 
It  reduces  for  a  ■  0  .  Laplace's  equation  is  known  to  have  5 
Independent  solutions  of  second  degree: 

22  -  (X2  +  y  2)/2  , 

xz,  yz,  ( 2- 17 ) 

(x2  -  y2  )/2,  xy  ; 

cf.  (Heiskanen  and  Moritz,  1967,  p.  61),  constant  factors  do 
not  matter.  Similarly,  Poincare's  equation  has  the  solutions 

z2  -  j  (1  -  4n2/o2  )(x2  +  y2  )  , 

xz,  yz,  (2-18) 

7  ( x2  -  y2  ) ,  xy  » 


as  we  see  by  substituting  Into  (2-16).  To  get  an  appropriate  time 
dependence,  these  solutions  must  be  multiplied  by  cosot  or 
sinot  . 

We 'shall  take  the  solution 

<j»  «  B(xz cosot + yz sinot)  ,  (2-19) 


which  is  a  linear  combination  of  the  solutions  xz  and  yz  . 
Using  (2-4)  and  (2-19)  in  (2-12)  we  get 
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(e2  -  4n2)u  ■  B( 1  -  2n/o }z  coset  , 

(az  -  4n2)v  ■  B(1  -  2fl/o)z slnot  , 
o2 w  ■  (B  -  2oQe)(x  cosot  *y  slnot)  . 

By  moans  of  the  substitution 

®  *  a\a  +  2fl) 
this  becomes 

u  «  8 z  cosot  , 
v  «  b  z  si  not  , 

w  ■  B  -  (x  cosot  +  y  slnot)  . 


(2-20) 


(2-21) 


For  e  «  0  ,  this  reduces  to  (1-16),  as  It  should. 

Actual  state  and  reference  state.  We  shall  Introduce  a 
"reference  state"  from  which  we  count  deformations.  It  corresponds 
to  the  absence  of  external  forces  or  deformations  ("undeformed 
state").  We  shall  denote  quantities  belonging  to  the  reference 
state  by  the  subscript  0  . 

Thus  we  may  split  up  the  total  potential  V  acting  on  a 
point  of  the  deformed  liquid  body: 


V  *  V0  +  Ve  +  Vi  . 


(2-22) 


V0  is  the  attractional  potential  of  the  undeformcd  body,  V# 
is  the  potential  of  external  forces  (tidal  potential),  and  Vi 
is  the  change  of  attractional  potential  of  the  body  due  to  the 
change  of  its  shape  caused  by  tidal  deformations. 

The  pressure  p  is  split  up  similarly  : 


P  *  Po  +  Pi  . 


(2-23) 


where  p  is  the  pressure  in  the  deformed  body*  and  po  is  the 
pressure  in  the  undeformed  body*,  pi  is  a  small  correction.  The 
constant  density  p  is  assumed  the  same  in  the  reference  state 
as  in  the  actual  state. 

Introduce  the  gravity  potential  Wo  in  the  undeformed  body 
as  the  sum  of  gravitational  potential  V0  and  centrifugal  poten¬ 
tial  e  : 


W0  -  V0  +  *  .  (2-24) 

Let  us  further  assume  that  the  reference  state  corresponds  to 
hydrostatic  equilibrium,  in  which  the  particles  are  at  rest.  Rota 
tion  is  about  the  z~axis,  so  that  e0  is  zero.  Then 

ju0  «  0 

since  there  is  no  particle  motion, 

♦  o  »  0 


by  (2-4)  since  c0 


0  ,  and  (1-38)  gives 


0  *  -gradij/Q 


(2-25) 


where 


-<l>o  *  V0  +  *  -  p-ipo  =  W0  -  p-ipo  (2-26) 

by  (2-2)  for  the  reference  state  (note  that  the  present  i|>0  is 
not  quite  the  same  as  (1-36)). 

The  substitution  of  (2-26)  into  (2-25)  gives 

grad  W0  *  £  grad  p0  ,  (2-27) 

which  is  the  well-known  condition  of  hydrostatic  equilibrium: 
cf.  (Jeffreys,  1962,  sec.  4.03).  The  physical  meaning  of  this 
equation  can  be  seen  by  writing  it  as 

dW0  =  £dp0  .  (2-28) 

Thus  d  po  =  0  implies  dW  *  0  and  vice  versa  :  the  surfaces  of 
constant  potential  coincide  with  the  surfaces  of  constant  pressure. 

Consider  now  such  a  surface  of  constant  reference  potential 
E o  ,  which  is  deformed  into  a  surface  z  (Fig.  2.1).  A  point  B 
of  the  actual  state  is  separated  from  the  corresponding  reference 
point  by  the  deformation  vector  u  . 
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n 


deformed  state 


reference  state 


Fig.  2.1.  Deformed  state  and  reference  state. 


Consider  now  any  quantity  f  of  the  actual  state,  and 
its  corresponding  quantity  f0  for  the  reference  state.  Differ¬ 
ences  of  the  form 

fA  -  f0A  or  fB  -  f0B  (2-29) 

in  which  f  and  f0  refer  to  the  same  point,  are  called  Eulerian 
di fferences ,  and  differences  of  the  form 

f8  -  foA  •  (2-30) 

in  which  f  refers  to  the  "actual  point"  B  and  f0  to  its  "re¬ 
ference  point"  A  ,  are  called  Lagrangian  differences.  This  is  an 
established  terminology  in  continuum  mechanics. 

To  the  geodesist,  an  analogy  to  the  gravity  field  may  be 
helpful.  The  actual  gravity  field  is  split  up  into  a  normal  field 
and  an  anomalous  field  (Heiskanen  and  Moritz,  1967,  sec.  2-13).  The 
surface  z  may  be  compared  to  the  geoid,  and  To  the  ellipsoid. 


ZSSSBS ElSSliS 
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The  anomalous  potential  T  ,  for  which  actual  potential  W  and 
normal  potential  U  refer  to  the  same  point,  is  an  "Eulerian" 
quantity,  whereas  the  gravity  anomaly  Ag  ,  for  which  actual 
gravity  refers  to  a  geodial  point  P  and  normal  gravity  y  to 
the  corresponding  ellipsoidal  point  Q  ,  is  a  "Lagrangian"  quant- 
i  ty . 

We  shall  now  make  the  definition  of  Vi  in  (2-22)  more 
precise  by  understanding  it  as  an  Eulerian  quantity  :  V  and 
the  reference  potential  V0  refer  to  the  same  point  : 

VA  *  V0A  +  VeA  +  V1A  .  (2-31) 

On  the  other  hand,  the  anomalous  pressure  pi  in  (2-23)  is  defined 
to  be  a  Lagrangi an  increment  : 


Pi  !  Pb  *  Poa  •  (2-32) 

Let  us  denote  by  n  the  difference 

"  *  “b  •  WA  • 

By  Taylor's  theorem  it  is,  to  the  linear  approximation  we  are  con¬ 
sistently  using, 

n  =  u  •  gradW  =  u  •  grad  W0  .  (2-33) 


Then 
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POB  *  POA  +  Pn  • 

Finally, 

Pi  3  Pb  "  Poa  *  Pb  *  Pob  +  p n  , 
so  that 

Pi*6p  +  p;i  ( 2-34 ) 

links  the  Lagrangian  increment  pi  to  the  Eulerian  increment 

«P  3  P  -  Po  .  (2-35) 

p  and  po  referring  to  the  same  point  B  .  (The  geodetic  reader, 
will  note  the  analogy  between  pi  and  the  gravity  anomaly  a g  ,  and 
between  <5p  and  the  gravity  disturbance  <$g  I) 

This  simple  but  important  relations  enable  us  to  transform 
expression  (2-2)  for  our  basic  function  : 

-♦  3  V  +  *  -  P-1p  +  (1  +  2.  ♦) 

*  V0  +  V e  +  Vi  +  *  -  p“1Po  -  P " 1  <5 P  +  (1+^4) 

3  W0  4  Ve  ♦  Vi  +  (l  +  £*)  -  „  .  (2-36) 

Here  we  have  used  (2-24),  (2-34),  and  (2-35).  The  integration  of 
(2-28)  gives 

Wo  *  ,  (2-37) 
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arart  from  an  integration  constant  which  we  can  put  equal  to  zero. 
The  use  of  (2-37)  in  '2-36)  gives  finally 

-*  *  V  +  n  -  ,  (2-38) 

if  with  Molodensky  we  put  * 

V  *  Ve  +  V !  +  (1  +£♦)  (2-39) 

(this  is  not  the  same  V  as  in  (2-22)  above!). 

External  and  induced  potential.  We  assume  an  external  (tidal) 
potential  of  the  form 

Ve(x,y,z)  *  k (xz  cosot  ♦  yz  sinot)  ,  (2-40) 

k  being  a  constant.  This  corresponds  to  a  tesseral  tidal  poten¬ 
tial  of  degree  2  and  order  1  ;  cf.  (Moritz,  1980b,  p.  30).  In 
fact,  transforming  to  spherical  coordinates  r,6,x  by 

x  *  r  si necos a  , 

y  *  r  sinesinx  ,  (2-41) 

z  =  r cosQ  , 

we  get 

ve('r'»9»x)  3  <  r2  s  i  necose  ( coscr  tcos  \  +  sinotsinx) 

*  < r2s i necosecos (ot  -  x )  . 
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With  the  abbreviation 

Si(e,x,t)  ■  sinecosacos (at  -  x)  (2-42) 

this  becomes 

V#  *  <r2Si  .  (2-43) 

We  readily  see  from  (2-4)  that  *  is  also  proportional  to  r2Sj 
Poisson's  equation  (Heiskanen  and  Moritz,  1967,  p.  5)  gives 

AVq  *  -4irGw  , 

where  a  denotes  Laplace's  operator  and  G  the  gravitational 
constant.  After  deformation,  the  density  p  remains  the  same. 
Thus,  also 

A(V0  ♦  VO  «  -4irGo  , 
whence  on  subtraction 

4VX  -  0  ,  '  (2-44) 

which  is  Laplace's  equation  for  harmonic  functions.  Hence  Vi  i 
harmonic  and  may  be  taken  to  be  proportional  to 

r2Si  ■  xz cosot  +  yz sinot  ,  (2-45) 


which  clearly  is  harmonic. 
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Hence  V  as  given  by  (2-39)  will  be  proportional  to  (2-45). 
The  expression  (2-19)  for  <i>  also  has  this  fcrm,  so  that  we  may 
put  with  Molodensky 

V  +  <|>  *  A(xz  cosot  +  yz  si  not)  .  (2-46) 

It  will  be  seen  that 

n  *  -C(xz  cosot  ♦  yz  sin  at)  (2-47) 

is  consistent  with  these  assumptions.  Inserting  (2-46)  and  (2-47) 
into  (2-38)  we  find  that  the  anomalous  pressure 

pi  ■  p  (A  -  C)  (xz  cosot  +  yz  sinot)  ,  (2-48) 

is  also  proportional  to  (2-45). 

A  boundary  condition.  The  normal  displacement  q  at  the  core 
mantle  boundary,  which  is  assumed  to  be  an  equi potenti al  ellipsoid, 
may  be  obtained  in  two  ways.  The  first  uses  the  fact  that  the  bound¬ 
ary  surface  is  an  equi potenti al  surface,  the  second,  that  it  is  an 
el  1 i psoi d. 

First,  let  Eo  and  E  in  Fig.  2.1  denote  this  boundary  sur¬ 
face  in  its  undeformed  and  deformed  position.  Then  the  unit  normal 
vector  £  is  given  by 


1  ■  -  ■  -S-VadU0  . 


(2-49) 


If 


g  -  |  grad  Wo  I 


denotes  gravity  (it  is  well  known  that  the  gravity  vector  grad  W0 
is  normal  to  the  equipotential  surface  £o  )•  Thus,  the  dot  deno¬ 
ting  the  inner  product, 


q  .  u  •  n  «  -g_1  u  *  grad  Wo  , 


(2-50) 


and  by  (2-33), 


q  a  -g"1 n 


(2-51) 


On  the  other  hand,  the  boundary  of  the  core  is  an  ellipsoid 
of  revolution 


F(x,y,z)  « 


x2  +  y2  .  z2 


+  -  1  •  0  , 


(2-52) 


a  and  c  denoting  its  semi-axes.  The  normal  vector  is  from  geome¬ 
try  given  by 


a-tWi 


(2-53) 


With 


.5  a2  -  c2 
e‘  ■  — $ —  , 

a4 
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e  denoting  the  excentri cl ty ,  we  have 

(1  -  *2)x 

grad  F  ■  X  (l-*2)y  .  (2-54) 

c 

z 

. 

To  a  sufficient  approximation, 

[grad  M2*  i  (xJ+y2  +  z2)  •  -y  4  Zi  (2-55) 

cH  c2  b2 

since  the  ellipsoid  is  very  nearly  a  sphere  of  radius  b  =  3470  km. 
Hence, 

'  (1  -  e2)x 

n  ■  e-  (l-e2)y  ■  (2-56) 

z 

Thus , 

q  *  u  •  n 

gives 

bq  ■  ( 1  -  e2 ) ( xu  +  y v )  +  zw  ,  (2-57) 


and  using  (2-21)  we  find 
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bq  »  ( 1  -  e2  )S  ( xz  cosat  +  yz  si  not  )  + 


g  -  12.  cl  (xz  coso  t  +  yz  si  no  t)  . 


Putting 


o  +  a  n  r. 

—  8  -  -  e 

or  a 


(2-38) 


this  becomes 


q  ■  b"1  (v  -  e2B)(xz  cosot  +  yz  sinot)  . 


(2-59) 


For  the  principal  lunisolar  tides,  o  •  -ft  (Moritz,  1980b, pp.  121- 
122),  so  that  v  and  v  -  e28  are  small  even  if  8  is  large,  which 
can  be  the  case  as  we  shall  see  in  sec. 10.  This  shows  that  the  appro¬ 
ximation  (2-55)  is  in  fact  justified. 

The  comparison  between  (2-51)  and  (2-59)  gives 


n  *  -gb* 1  (\>  -  e28 )  (xz  cosot  +  yz  si  not)  . 


(2-60) 


The  gravity  g  at  a  point  of  radius  vector  r  at  the  interior  of 
a  homogeneous  sphere  is  well  known  to  be 


g  »  pr 


(2-61) 


(Kellogg,  1929,  p.  19).  Putting  r  *  b  and  substituting  into 
(2-60)  we  get 


,s  i...^grai 


—  tn 


n  ■  -  p  (v  -  e2S)(xz  co sot +  yz  si  not)  . 


(2-62) 


The  comparison  with  (2-47)  gives 

C  -  P(v  -  e28)  .  (2-63) 

3 

If  we  Introduce  a  new  constant  y  by 

c  -  -  ili  Pe2y  ,  (2-64) 

then  (2-63)  gives 

(y  - B)e2  +  v  «  0  .  (2-65) 

which  can  be  considered  a  boundary  condition  expressing  y  in  terms 
of  8  and  e  in  view  of  (2-58). 

The  constants  A  and  a  .  By  (2-46)  we  have 
V  +  ij»  »  A(xz  cosot  +  yz  sinot)  . 

Instead  of  A  ,  let  us  introduce  another  constant  a  ,  which  is 
related  to  A  in  the  same  way  (2-64)  as  C  is  to  y  • 
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This  parameter  o  was  used  by  Molodensky  (1961)  who  has  put 


V  +  ♦  ■  o  ♦  , 

*  «  K(r)r“2(xz  cosot  ♦  yz  slnot)  . 


(2-67) 


The  function  K(r)  satisfies  a  differential  equation  which  for  a 
homogeneous  core  becomes 


r  +  -  K'  -  JL  k  >»  o 

r  rz 


(2-68) 


A  solution  is 

K(r)  -  Dr2  .  (2-69) 

Putting  the  constant  0  ■  Aa~l  ,  we  get  (2-46). 
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3.  Toroidal  and  Spheroidal  Modes 


The  modern  treatment  of  core  motion  Is  In  terms  of  toroidal 
and  spheroidal  oscillations  (Shen  and  Mansinha,  1976;  Smith*  1977; 
Wahr*  1979).  We  shall  relate  Molodensky's  solution  to  this  approach 
Toroidal  and  spheroidal  modes  are  well  known  from  the  theory 
of  free  oscillations  of  the  earth;  cf-.  (Kaula*  1968,  sec.  2.2). 

They  were  originally  considered  for  an  elastic  body  but  can  be 
applied  to  a  liquid  body  as  well,  since  the  latter  may  be  regarded 
as  a  limiting  state  of  an  elastic  body;  cf.  sec.  5. 

The  displacement  vector  u  Is  expressed  as 


“  ■  n  <sr  ♦  i”)  •  (3-D 

nm 

S“  and  T®  denoting,  respectively,  spheroidal  and  toroidal  modes 
of  degree  and  n  and  order  m  . 

In  spherical  coordinates  r,e,x,  these  vectors  S®  and  T® 
have  the  components 


(S“)r  *  H®(r)Pnm(cose)cos(ot  -  mx)  , 

(s“)6  ■  L"(r)  ±  P„»(cose)eos(ot  -  mx }  .  (3-2) 

(i”)x  ■  ■>»-“(»->  r(W  Pno(cos«)stn(»t-»ix)  ; 


(£),  • 0  • 

(£),  •  -i»T“(r)  Tnis-  P^lcosejtosf’t  -  *M  , 
(Six  *  -S<r>  TC  Pnm(«s«)sU(at-«»)  . 


(3-3) 


We  have  followed  the  real  notation  of  Shen  and  Mansinha  (1976); 
other  authors  use  a  complex  notation  with 

eiot  ■  cosot  1  si  not  .  (3-4) 

The  quantities  and  T*  are  functions  of  the  radius  vector 

r  ,  and  pnmUos«)  ere  the  usual  Legendre  functions;  cf.  Heiskanen 
and  Moritz,  1967  sec.  1-11) . 

We  consider  only  and  T J .  With 

Pu(cose)  «  sine  ,  P2l(cose)  «  3sinecose 

we  get  from  (3-2)  and  (3-3) 

ui)r  •  o . 

(l})e  «  -Tcos(ot-  x)  ,  (3-5) 

(t{)  ■  -Tcosesi n(ot  -  x  )  ; 

(S 2)  •  3Hcosesinecos(ot  -  x)  , 

(3-6) 


1  (3-6) 

(S2)0  *  3L(cos2e  -  sin2e  )cos(o  t  -  x  )  , 

(S2)x  *  3Lcose  si  n(o  t  -  x  )  , 
where  we  have  abbreviated 


H 2 ( ^ )  *  H(r)  =  H  , 

L2 ( r )  =  L(r)  =  L  ,  (3-7) 

T}(r)  =  T(r)  *  T  . 


The  components  of  the  vector  jj  in  spherical  coodinates  are 
transformed  into  its  components  Uj,  u2  ,  u3  in  rectangular  coordi¬ 
nates  x.y.z  by  a  rotation  : 


u 

u 

u 


1 

2 
3 


!  sinecosx 
si nesi nx 
cose 


i. 


m 

-  *i 

cosecosx  -sinx 

u  ! 

r 

cosesinx  cosx 

ue 

-sine  0 

_ux_ 

(3-8) 


cf.  (Heiskanen  and  Moritz,  1967),  p.  206  ,  equs .  (5-54a,b,c),  putting 
e  =  90°  -  <}, 

If  we  retain  only  toroidal  oscillations  (1,1)  and  spheroidal 
oscillations  (2,1),  then  (3-1)  reduces  to 


u 


(3-9) 


and  both  vectors  on  the  right-hand  side  are  transformed  by  (3-8), 
so  that  we  get  rectangular  components  from  spherical  ones. 


using 


By  a  straightforward  though  somewhat  laborious  computation, 
(2-41) ,  we  obtain 

(tJ)  1  *  -Tr-1  z  coso t  , 

(lj)a  c  "Tr"1  z  sinot  ,  (3-10) 

(Tj)3  «  Tr"1  (x  coso  t  +  .y  si  no  t)  ; 

(S^)  i  *  3(H  -  2L)r“3xz(x  cosot  +  y  sinot)  +  3Lr-1z  cosot  , 

(S1)  n  *  3(H  -  2L)r"3yz(x  cosot  +  y  sinot)  +  Sir”1  z  s i  no  t  , 

■**  4 

(3-11) 

( S_2 )  3  *  (  3 ( H  -  2t)r“3z2  +  3Lr-1]  (x  cosot +y  sinot) 


The  solution  (2-21)  may  be  simplified  by  means  of  (2-58)  : 


U  »  BZ COSot  , 
v  =  Bz  si  not  , 

w  *  (v  -  6  )(x  cosot  +  y  sinot)  . 

This  may  be  written  as  follow;: 
u  ■  +  jj"  ; 

U '  ■  ( B  -  Y  v )  z  cosot  , 

v'  a(B-Yv)zs‘n°^  » 

w1  *  -  ( 8  -  ^  v  )  ( x  cosot  +  y  si  not )  ; 


(3-12) 


(3-13) 


(3-14) 
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U "  *  ^  v  Z  COSO  t  , 

v"*-£vzsinot  * 

w1'  *  i  v(x  cosot  +  y  si  not)  . 

The  comparison  with  (3-10)  and  (3-11)  shows  that 

u'  ■  I}  ,  u"  ■ 

with  H  ■  2L  and 

H  ■  ?  vr  , 

L  »  ^  vr  , 

T  ■  -  (B  -  \  v)r  , 
where,  by  (2-58) 


(3-15) 


(3-16) 


(3-17) 


(3-18) 


This  shows  that  Molodensky's  basic  solution  (2-21)  may  be  inter¬ 
preted  as  the  resultant  of  a  toroidal  and  a  spheroidal  oscillation, 
the  functions  (3-7)  being  simply  proportional  to  r  by  (3-17). 

Another  insight  is  provided  by  considering  the  matrix  of  defor¬ 
mation  gradients 
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U 


Here  we  have  put 

t 


*  x  , 


X2  -  y  . 


z 


and  denoted  partial  derivatives  by 


_ 


3  X  ’ 


The  matrix  U  can  be  split  up  into  a  symmetric  part  l 
antisymmetric  part  R  : 

U  •  E  +  R  , 


where 


—  * 


■ij 


1 

7 


3Ui  3U^ 


3X. 


3  X . 


j 

is  the  strain  tensor  (see  also  sec.  5)  and 


tru]  • 


r . 


1 


3  "  7 


3  U 


3  X 


3 


3  U 
3  X 


is  the  (infinitesimal)  rotation  tensor. 


(3-19) 


(3-20) 


and  an 

(3-21) 


(3-22) 


(3-23) 
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For  the  solution  (3-12)  we  have 

0  0  Bcosot 

0  0  Bsinot  ,  (3-24) 

( v  -  B  )  coso  t  (v-B)slnot  0 


0 

0 

1 

7 

VCOSO  t 

0 

0 

1 

7 

\> si  not  , 

(3-25) 

vcosot 

^  v  si  not 

0 

0  (B  -  y  v  )  CO  So  t 

0  (B--ij-v)sinot 

( B  -  v  )  s  1  no  t  0 

(3-26) 

The  comparison  with  (3-14),  (3-15),  and  (3-16)  shows  that 

u1  -  l\  *  R x  ,  u"  =  S‘  -  E x  .  (3-27) 


The  symmetric  strain  tensor  E  ,  which  characteri zes  an  affine  de¬ 
formation  of  the  volume  element,  is  thus  seen  to  be  related  to  the 
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spheroidal  mode  (2,1),  and  the  antisymmetric  rotation  tensor  R 
is  related  to  the  toroidal  mode  (1,1),  which  thus  represents  an 
(infinitesimal)  rotation. 

The  amplitude  of  this  rotation,  as  a  function  of  r  ,  is  re¬ 
presented  by  T(r)  ,  whereas  H(r)  characteri zes  radial  deforma¬ 
tion  and  l(r)  ,  tangential  deformation  (displacement  normal  to 
the  radius  vector).  This  is  immediately  evident  from  (3-2). 
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PoincarS  (1910,  sec.  II)  expressed  the  coordinates  x,y,z 
of  a  fluid  particle  as  a  linear  function  of  its  coordinates  x0 , 
y0,z0  in  a  standard  position: 

x  ■  #1 1 ( t)x0  +  a12(t)y0  +  a13(t)z0  , 

y  *  a2i(t)x0  ♦  a22(t)y0  +  a23(t)z0  ,  (4-1) 

z  *  a31(t)x0  +  a32 (t)y0  +  a33(t)z0  , 

or  briefly 

xi  “  aij(t)xoi  »  (4‘2) 

using  the  notation  (3-20),  the  summation  convention  implying  sum¬ 
mation  over  a  subscript  occurring  twice  (in  our  case,  j  ). 

Jeffreys  (1949)  used  this  method,  employing  the  "nutation 
frame"  x°x2x3  (sec.  1).  Thus  we  should  have  to  write  x°  in¬ 
stead  of  xL  ,  and  Xoi  instead  of  x0i  ,  but  we  shall  take 
the  superscript  0  for  granted.  Thus,  Xoi  are  the  coordinates 
of  a  liquid  particle  in  an  "undisturbed  reference  state",  in  which 
the  liquid  is  at  rest  in  the  nutation  frame  (rotating  uniformly 
with  respect  to  inertial  space). 

Then  the  displacement 
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will  be  small,  and  (4-2)  gives 


xi  “ 


Oi 


ai  j  X0  j  '  X0i 


or 


01  X 

ij  03 


(4-4) 


wi  th 


a 


ij 


(4-5) 


the  Kronecker  delta  6^  denoting  the  elements  of  the  unit  matrix 
as  usual  . 

The  solution  (2-21)  is  clearly  of  this  linear  form.  Generally, 
the  linear  representation  (4-4)  is  seen  to  be  related  to  a  quadratic 
pertu:‘  ng  potential  ,  a  quadratic  pressure  pi  ,  and  a  quad¬ 
ratic  function  <ji  as  (2-12)  shows.  In  terms  of  spherical  harmonics, 
it  corresponds  to  perturbati ons  which  are  harmonics  of  che  second 
degree . 

The  transition  from  the  standard  position  to  the  actual 

positon  x  may  be  described  in  the  following  way  (Fig.  4.1).  First, 
there  is  a  dilatation  which  transforms  the  elliptical  core  into  a 
sphere*  it  '  represented  by  the  matrix 


1  0  0 


Dilatation  Q 


Compression  Q'1 


FIGURE  4. 


Method  of 
Jeffreys 


Rotation  R ' 
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where  a  and  c  denote  the  semiaxes  of  the  core,  considered  as 
an  ellipsoid  of  revolution.  The  next  step  is  a  small  rotation  by 
the  matrix 


R  *  I  + 


0  -n 
n  0 


1 

m 


■  I  +  F 


!  -1  -m 
c 


0 


(4-7) 


which  differs  from  the  identity  2.  the  sraa11  antisymmetric 
matrix  F  .  The  compression  by  the  matrix  D"1  transforms  the 
sphere  back  into  our  ellipsoid,  and  a  final  rotation 


R'  -  I  + 


0  -n'  r 
n‘  0  m' 
-T  -m'  0 


I  +  F' 


(4-8) 


gives  the  actual  position  (m'  has  no  relation  to  m  ,  etc.). 
The  composition  of  these  linear  transformations  gives 

*.  *  R '  D’ 1  R  0  x0  ■  (2.+  4>)*o  *  £0  +  ££0  »  ( 4-9 ) 

where  the  small  matrix  a  is  expressed  by 
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A  «  (I  ♦  F' )D‘1(A  +  F)D  -  l  .  (4-10) 

Retaining  only  terms  linear  in  F  and  F'  we  get 

A  *  *  F ’  +  D" 1  £_£  ;  (4-11) 

in  fact,  (4-3)  and  (4-9)  give 

u  -  a  * o  •  (4-12) 

so  that  the  comparison  with  (4-4)  shows  that  the  matrix  A  has 
the  elements 

Working  out  the  matrix  multiplications  in  (4-11)  we  get,  ex 
eluding  quadratic  and  higher  terms, 


“11  ■ 

“22  *  ®33  *  0  » 

a  1  2  * 

- (n  ♦  n 1 )  , 

a2 1  *  n 

+  n'  , 

“13  = 

f1  +  • 

i 

“31  *  - 

[f1  +1‘1 

» 

(4-13) 

“2  3  * 

|m  +  m1  , 

“32  »  ~ 

I™ * "'■)  ■ 

Note 

that  the  vectors 

x,x  ,u 

refer  to  the 

nutation  frame 

Equations  of  motion.  Let  us  use  the  equation  of  motion 
(1-13)  for  the  nutation  frame: 
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««  » 
u  -  2Qv 

v  +  2Qu 

•« 

w 


li 

ax  * 
ii 

a  y  » 

a+ 

n  * 


(4-14) 


ai.d  substitute  Jeffreys'  representati on 


ui  "  aijxoj  » 

in  which,  of  course,  u^  «  u  ,  U2  »  v 
i nstance , 


(4-15) 

U3  »  w  .  We  have,  for 


u  a  “ljXOj  *  auxo  +  “l2*0  +  au2o  •  (4"l6> 

and  the  differentiation  gives 

since  the  reference  position  x0^  does  not  change  with  time. 
Thus  the  system  (4-14)  becomes: 


°ljXOj 
°  2  j  X0 j 
°  3  j  X0 j 


-  2Pa 
2no 


2jX0j 

ljXOj 


(4-18) 


ik 


Since  a 


and  its  derivatives  are  small,  we  may  put 
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Xoi  4  *1  *  X  , 

x02  4  x2  -  y  ,  (4-19) 

X0  3  *  X3  -  Z 

In  these  equations.  Let  us  now  differentiate  the  first  equation 
of  (4-18)  with  respect  to  z  and  the  third,  with  respect  to4  x. 
In  view  of  (4-19)  we  get 


a  1  3  “  2&a2  3 


3X3Z 


(4-20) 


We  subtract  both  equations,  obtaining 

»*  ••  * 

“13  ’  °31  "  2Qa23  a  0  .  (4-21) 

From  the  second  and  third  equation  of  (4-18)  we  get  in  the  same 
way 


••  ••  • 

a23  *  a32  +  Zftaj  3  *  0  . 


(4-22) 


(4-21)  and  (4-22)  represent  equations  of  motion  in  terms  of  matrix 
elements  . 

In  the  remaining  part  of  this  section  we  shall  compare  results 
of  the  methods  of  Poincar§,  Jeffreys,  and  Molodensky.  To  make  essen¬ 
tial  relations  transparent,  we  shall  limit  ourselves  to  Poincarfe's 
model :  the  mantle  is  riaid,  the  earth's  surface  and  the  core-mantle 
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boundary  are  concentric  and  coaxial  ellipsoids  of  revolution, 
and  the  core  consists  of  a  homogeneous  and  incompressible  liquid 
This  is  the  model  considered  In  (Morit2,  1980b,  secs.  12  and  13) 
Relations  between  body  frame  and  nutation  frame.  The  co¬ 

ordinates  of  any  point  in  the  body  frame,  forming  the  vector  x 
and  of  the  same  point  in  the  nutation  frame,  forming  the  vector 
x°  ,  are  related  by  (1-17)  and  (1-18)  : 

x  »  (I  +  e)x°  ,  (4-23) 


0  « 


0 

-e3 

e2 


®  3 
0 

-6  x 


-e2 

«i 

0 


(4-24) 


Thus  1+0  rotates  from  the  nutation  frame  to  the  body  frame, 
whereas  the  matrix  (4-8)  rotates  the  body  frame  represented  by 
the  axes  of  the  ellipsoid)  to  the  inertial  frame.  Thus 


F' 


* 


-0 


(4-25) 


and  the  comparison  between  (4-8)  and  (4-24)  gives 
0 1  ■  -in '  ,  02  ■  1  1 


(4-26) 


We  shall  assume 
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9  3  *  n '  *  0 


(4-27) 


for  polar  motion  (Moritz,  1980b,  p.  85). 
We  introduce  the  complex  quantity 


4  w  =*  ej  +  i  e2  (4-28) 

(ibid. ,  p.  80),  and  define  similarly 

X  '  =  1 '  +  im '  .  (4-29) 

Then  (4-26)  is  equivalent  to  the  simple  complex  equation 

w  *  i X  *  .  (4-30) 

A  geometrical  interpretation  of  the  complex  number  x'  is 
obtained  by  considering  the  nutation  np  of  the  figure  axis 
(our  x3  axis).  We  have  (ibid.,  p.  128)  : 


np  *  -iw  , 
and  by  (4-30) , 


n  •  -  X '  .  (4-31) 

F 

Thus  x'  represents  simply  the  nutation  of  the  figure  axis. 

Let  us  now  consider  the  direction  of  the  rotation  axis  in 
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the  body  frame,  which  represents  nothing  else  than  polar  motion 
PR  ■  m!  +  im2  ;  (4-32) 

cf  .  (ibid.,  pp.  93  and  129).  By  (l-19a,b)  we  have 


u  =  dig  +  6w  * 


in  view  of  (l-19c),  where 


(4-33) 


mi  *  Qmi  ,  w2  a  flm2  .  (4-34) 

In  agreement  with  (ibid.,  p.  107)  we  put 

u  3  +  iu2  ,  (4-35) 

so  that 

u  »  Q(mi  +  im2)  =  npR  .  (4-36) 

The  relation  between  u  and  w  is 


u  *  w  +  inw 


(4-37) 


(ibid. ,  p.  115),  and  for  w  proportional  to  e 


as  usual  , 


w  *  woe 


lot 


(w0  =  const.) 


this  becomes 


u  *  i  (o  +  fl  )w  . 


Together  with  (4-36)  and  (4-30)  this  gives 


.  o  +  £2 

Pd  =  1  -o—  W 


q  +  fl  ,  • 


, 


in  agreement  with  ibid.,  eq ,  (13-38). 
By  (1-28)  we  have 


p  *  e  (cosot +  isinot)  =  ee 


io  t 


so  that 


ee 


lot 


( 1  ’  ♦  1  m ' ) 


(4-3 


(4-3 


(4-4 


(4-4 


(4-4 


represents  a  relation  between  Molodensky's  parameter  e  and 
Jeffreys'  parameters  1*  and  m’  . 

Relation  between  Jeffreys'  and  Poincarfe's  methods.  Le 
us  substitute  (4-13)  into  (4-21)  and  (4-22),  taking  account  of 


the  fact  that  a  and  c  are  constant  for  a  rigid  mantle.  The 
result  is 


|  +  i]  1  +  21'  -  2n  i  m  -  2flm'  *  0 

lc  aJ  c 

||  +  |J  m  +  2m'  +  2n  |  1  +  2 ft  1  '  -  0  . 


( 4-43 ) 


These  are  the  equations  of  motion  for  Jeffreys'  parameters. 

It  is  easy  to  bring  them  to  a  form  familiar  from  (Moritz, 

1980b,  sec.  12).  By  (ibid. ,  p.  103)  we  have  for  rotational  symmetry 

(a  =  b)  ' 

Ac  “  i  Mc<®2  +  c2)  . 

Cc  -  |  Mca2  ,  (4-44) 

F  »  Mcac  , 

Mc  denoting  the  mass  of  the  core  and  Ac  and  Cc  its  principal 
moments  of  inertia;  F  is  an  auxiliary  quantity.  Then,  on  multi¬ 
plying  them  by  Mcac/5  ,  eqs.  (4-43)  may  be  written: 

A  1  +  FI'  -  C  Jim  -  Film'  ■  0  , 

c  c  (4-45) 

A  m  +  Fm'  +  C  fll  +  Fnl '  «  0  . 

0  c 

This  system  has  been  used  by  Jeffreys  (1949,  p.  677,  eqs.  (13)  and 
(14)). 


Hcaatt 
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It  is  similar,  though  not  identical,  to  Poincar&'s;  equations 
(12-19)  in  (Moritz,  1980b,  p . 1 07 ) .  It  is,  however,  identical  to  a 
related  equation,  namely  to  the  second  equation  of  the  system  (12-46), 
ibid. ,  p.  115. 

To  see  this,  let  us  go  to  complex  notation,  using  (4-29) 
and,  similarly, 

X  *  1  +  im  .  (4-46) 

Then  the  system  (4-45)  is  equivalent  to  the  single  complex  equation 

A  X  +  Fx'  +  iC  QX  +  iFflX'  »  0  .  (4-47) 

o  e  '  ' 

Now  we  use  (4-30)  and  put,  similarly, 

v  «  ix  .  (4-48) 

Then  (4-47)  becomes 

••  •  ■ 

F(w  +  inw)  +  A  v  +  iC jjv  =  0  ,  (4-49) 

c  c 

identical  to  the  second  equation  of  (12-46),  ibid. 

The  fact  that  we  have  arrived  at  this  equation  in  a  completely 
different  way,  serves  as  a  check  and  shows  a  relation  between  Jef¬ 
freys'  (1949)  method  and  "Poincarfe's  first  method"  (Poincarfe,  1910, 
sec.  II),  of  which  Jeffreys'  method  is  a  modification  and  which  thus 
is  quite  similar,  except  that  Poincarfe  uses  al3,  a3<  ,  a23,  a32 
instead  of  1,  1',  m,  m'  ;  these  parameters  are  related  by  (4-13). 

The  equations  of  motion  for  Poincarfe's  parameters  are  (4-21)  and 
(4-22). 


. .  T. 3T&1* li I  .Mli  '  f  I. 


f  i  Sf!TOB»!W 


”*mm ppw 
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Relation  between  Molodensky's  and  Poincarfe's  methods. 

Using  the  complex  numbers  u  and  v  ,  given  by  (4-35)  and  (4-48), 
we  may  write  Poincarfe's  basic  equations  in  the  form  (Moritz,  1980b, 
eq.  (12-22);  we  take  L  *  0)  : 

Au  +  Fv  -  i(C-A)Qu  +  iFnv  *  0  ,  (4-50) 

•  * 

Fu  +  A  v  +  iC  Av  *  0  ,  (4-51) 

c  c 

A  and  C  denoting  the  principal  moments  of  inertia  of  the  whole 
earth . 

With 


u 


u0e 


tot 


v 


v0e 


iot 


(4-52) 


and 


(ibid, »  p.  109)  this  becomes 

[  Ao  -  (C-A)flluo  +  Ac(o  +  £l)vo  s  0  ,  (4-53) 

A  ouq  +  (A  o  +  C  JJ)vq  ■  0  ,  (4-54) 

c  c  c 

We  shall  now  show  that  PoincarS's  equation  (4-54)  is  a 
simple  consequence  of  Molodensky's  theory  applied  to  the  Poincarfe 
model  . 

We  use  Molodensky's  solution  in  the  form  (3-12)  : 
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Ui  »  8  X3  COSo  t 
U2  *  6X3Sinot 


(4-55) 


The  third  equation  we  shall  not  need.  These  equations  refer  to 
the  body  frame. 

These  equations  may  be  compared  to  (4-4).  There,  however, 
the  nutation  frame  is  used.  It  differs  from  the  body  frame  by  the 
final  rotation  (4-8).  If  this  rotation  were  not  performed,  we  should 
end  up  in  the  body  frame.  Since  this  is  desired,  we  just  do  not 
perform  the  final  rotation  or  simply  put  1  *  =  m '  =  n '  =  0  .  Then 
(4-4)  with  (4-13)  gives 


ui  *  f  lx3  » 
u2  »  I  mx3  ; 


(4-56) 


we  have  identified  X03  on  the  right-hand  side  with  X3  ,  which 
is  permissible  because  of  the  smallness  of  1  and  n  . 

The  comparison  between  (4-55)  and  (4-56)  gives  immediately 


1  *  8cosot 

Q 


m  »  £  esinot  , 

a 


(4-57) 


or,  with  (4-46 ) , 


(4-58) 


By  (4-48)  we  get 
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v 


Liot 


» 


and  by  (4-52) 


vo 


I 


(4-59) 


(4-60) 


which  is  a  relation  between  Poincar&'s  variable  v  and  Molodensky's 
parameter  8  . 

From  (4-36)  and  (4-41)  we  get  an  analogous  relation 

uo  *  Oe  .  (4-61) 

Thus  Molodensky's  parameters  e  and  6  are  seen  to  be  basically 
equivalent  to  u  and  v  . 

Let  us  now  consider  eq.  (2-59)  for  Molodensky's  theory. 

Since  the  core-mantle  boundary  is  assumed  to  be  rigid,  there  is  no 
normal  displacement: 

q  *  0  .  (4-62) 

Thus  (2-59)  gives  simply 

v  -  e 2  6  *  0  ;  (4-63) 

we  could  have  obtained  this  condition  also  from  (2-65).  It  will 
now  be  shown  that  (4-63)  is  equivalent  to  (4-54). 

First,  (2-58)  and  (4-63)  give: 


a  e 

o 


1 

7 


e2  8 
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or 


6+~(S-s)a^e26  . 
By  definition, 


and  the  comparison  of  eqs.  (12-2 
p.  109)  shows  that 


This  is  substituted  into  (4-64). 
lations  this  is  brought  into  the 

Ajoe  +  (V  +  Ccn  He  - e )  * 
Let  us  finally  bring  (4-60)  into 


(4-64) 

(4-65) 

)  and  (12-30)  of  (Moritz,  1980b, 

(4-66) 

By  elementary  algebraic  manipu- 
f  orm 

0  .  (4-67) 

a  different  form.  Using 

(4-68) 


we  get: 


-Vo  B  oP  -  £  o  e2p 


<*P  “ 


(by  (4-63) ) 


*  op  -  [  (o  +  Q)P  -  fie]  (by  (2-58)) 


which  gives  simply 


Vo  ■  o(P  -  c ) 


Now  (4-61)  and  (4-69)  are  substituted  in  (4-67)  and 
immediately  give  (4-54),  which  was  to  be  shown. 

In  sec.  10  we  shall  in  a  similar  way  derive  Poincare's 
equation  (4-53)  from  Molodensky's  theory,  completing  the  proof 
of  equivalence  of  both  theories  for  the  Poincare  model. 
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PART  B 


THEORIES  OF  MOLODENSKY  TYPE 


5 •  Equations  of  Elasticity 


The  strain  tensor.  We  briefly  summarize  basic  formulas 
from  elasticity  theory,  for  details  referring  the  reader  to 
any  textbook  on  theoretical  physics  or  to  special  treatises 
such  as  (Love,  1927)  or  (Green  and  Zerna,  1968). 

Let  x0  and  x  denote  the  position  vectors  of  a  point 
of  an  elastic  body  before  and  after  deformation  (Fig.  5.1). 
Then  the  difference 


FIGURE  5.1. 


Elastic  deformation. 
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u  -  £  *  x0  (5_1) 

is  the  displacement  vector,  or  deformation  vector.  Denoting  the 
components  of  jj  by  ui  as  usual ,  we  get  the  matrix  of  defor¬ 
mation  gradients 


U 


(5-2) 


in  view  of  the  smallness  of  ui  it  is  irrelevant,  to  the  linear 
approximation,  whether  we  differentiate  with  respect  to  x0i  or 
to  x  ^  . 

The  symmetric  part  of  U  , 


E  « 

1 

7 

(M 

*  UT' 

^  i 

|  *  Ceij] 

1 

1 

l-k 

+  ila1 

e  ij 

* 

7 

» 

.  9  *  i 

3Xi 

(5-3) 

(5-4) 


is  the  strain  tensor,  and  the  antisymmetric  part  of  U  , 


(5-5) 

(5-6) 


is  the  (infinitesimal)  rotation  tensor. 


All  this  is  fully  ana 
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logous  to  sec.  3. 

The  rotation  tensor  R  plays  rather  a  small  role  in  elasti¬ 
city.  0,'  basic  importance  is  the  strain  tensor  E  .  Denote  a 
line  element  before  and  after  deformation  by  dso  and  ds  ,  and 
a  volume  element,  by  dvo  and  dv  .  Then  the  change  of  ds  is 
given  by  , 

ds2  -  ds§  «  2dxTEdx  ,  (5-7) 

and  the  volume  dilatation  3  by 

©  «  dv-~-  d-vQ-  ■  trace  E  •  eu  +  e22  +  e33  » 

dv  - 


Hence  the  change  of  both  length  and  volume  can  be  expressed  in 
terms  of  the  tensor  l  . 

We  shall  now  again  use  the  Einstein  summation  convention: 
if  an  index  occurs  twice  in  a  product  or  a  similar  expression,  sum¬ 
mation  with  respect  to  that  index  is  automatically  implied.  For 
i nstance , 

*Jbi  ■  j/l*)  ■  Vi  +  +  V3  • 


(5-9) 
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and  (5-8)  reduces  to 


0  »  e 


kk 


(5-11) 


The  stress  tensor.  There  are  two  basic  tensors1  in  ela¬ 
sticity  theory.  The  first  is  the  strain  tensor  £  ,  the  second 
the  stress  tensor,  denoted  by 


P  -  tPij]  .  (5-12) 

This  tensor  is  characterized  by  the  fact  that  the  stress  vector 
£  on  any  surface  element  normal  to  a  unit  vector  n.  is  expres¬ 
sed  by 


£  3  ££ 


(5-13) 


or 


Pi  »  Pijnj  .  (5-14) 

The  tensor  P  is  symmetric: 

PT  *  P  .  Pji  ■  Ptj  •  (5-15) 


n  A  tensor  may  be  viewed  simply  as  a  matrix  having  physical 
meaning,  such  as  inertia  tensor,  strain  tensor,  or  stress 
tensor,  Mathemati cl ans ,  of  course,  characterize  a  tensor 
in  terms  of  its  behavior  with  respect  to  coordinate  trans¬ 
formations,  but  this  will  not  be  needed  here. 
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An  important  special  case  is  that  of  P.  being  diagonal 
and  proportional  to  the  unit  matrix  I_  : 

p  3  -Pi  Pi;j  -  -P^ij  *  (5-16) 

8  being  the  Kronecker  delta  already  used  in  sec.  4  :  <5^  =  1 

or  0  according  to  whether  j  *  i  or  j  t  i  ;  it  denotes  the 
elements  of  the  unit  matrix.  The  proportional i ty  factor  p  is  called 
hydrostatic  pressure  since  the  form  (5-16)  plays  a  special  role 
in  hydrodynamics. 

Strain  tensor  i  and  stress  tensor  P  are  linearly  rela¬ 
ted  by  Hooke's  law: 


PlJ  1  XekkSlj  *  Z"ei3 
or,  by  (5-11) 

P  =  AOj_  +  2uE  , 


(5-17) 


(5-18) 


where  A  and  u  are  Lamfe's  parameters.  In  terms  of  components 
we  have,  for  instance, 


Pn  =  XQ  +  2y  lit  * 

5U  +  ^2  1 
,  3X2  j 


P 1 2  =  V 


(5-19) 
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Equations  of  motion.  They  have  the  form 


p u.  *  p  f .  + 
i  i 


(5-20) 


where  p  is  the  density,  fi  is  the  body  force  per  unit  mass, 
and  3pij/3xj  represents  the  resultant  force  due  to  the  stresses 
p.  acting  on  a  volume  element.  In  vector  notation  we  may  write 
(5-20)  as 


pjj  =  pf  +  divP  .  (5-21) 

The  substitution  of  Hooke's  law  (5-17)  gives 

p-»i  -  + »*)'■—■  +  uAUi  +  pft  ,  (5-22) 

or 

p^  a  (x  +  p)grade  +  yau^  +  p_f  ,  (5-23) 

©  being  the  volume  dilatation  (5-11)  and 

32U.  3  2  U  .  3  2  U  .  3  2  U  . 

A  Ui  =  - —  =  - —  +  - -  +  - -  .  (5-24) 

3X  3X  3Xf  3xg  3  X^ 

denoting  Laplace's  operator  applied  to  the  displacement  component 

Ui  . 

Relation  to  hydrodynamics.  Hydrostatic  fluid  pressure  has 
the  form  (5-16).  The  comparison  of  (5-16)  and  (5-18)  shows  that, 
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for  this  case, 

w  =  0  ,  X 0  =  - p  ,  (5-25) 

Then  (5-23)  reduces  to 
p  u  *  -grad  p  +  p  f 


or 


jj  *  f  -  p_1grad  p  ,  (5-26) 

••  ** 

which  is  (1-1)  with  _u  =  and  f_  «  grad  V  . 

Thus  we  have  shown  the  basic  fact  that  the  equations  of 
hydrodynamics  may  formally  be  considered  a  special  case  of  the 
equations  of  elasticity  if  the  Lamfe  constants  are  replaced  by 
(5-25).  For  incompressible  fluids,  there  is  0=0  but  p  fini¬ 
te,  so  that  for  them  we  must  take 

u  =  0  and  X  -►  «®  ( 5-27 ) 

in  such  a  way  that 

lim(X0)  «  -p  (5-28) 


remains  finite. 
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6 .  Application  to  the  Earth 


In  order  to  apply  the  elastic  equations  of  motion  (5-20) 
or  (5-21)  to  the  rotating  earth,  we  must  transform  them  to  a 
system  xyz  rotating  with  the  constant  angular  velocity  vector 
w  .  In  this  system  we  have  as  in  sec.  1  : 

x  =  x0  +  u  , 

x,  *  u  .  (6-1) 

••  «« 

X  *  u  , 

denoting  an  undisturbed  equilibrium  state  corresponding 
to  the  uniform  rotation.  Then,  in  complete  analogy  to  (1-4), 
eq.  (5-21)  is  replaced  by 


••  • 

p  [  u^  +  2u  *  +  u  x  (u  xx,  )  ]  =  pgradV  +  div£  .  (6-2) 

The  comparison  of  the  left-hand  sides  of  (1-4)  and  (6-2)  shows 
that  we  have  used  (6-1)  and 

u  «  0  (6-3) 

for  uniform  rotation.  Note  that  this  means  that  the  coordi nate 
system  xyz  ,  not  the  earth,  rotates  uniformely:  small  non-uni 
formities  of  the  rotation  of  an  elastic  earth  are  included  in 
the  vector  u^  . 

Furthermore  we  have  put 
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f  *  grad  V  ,  (6-4) 

expressing  the  gravitational  force  as  the  gradient  of  the  gravi¬ 
tational  potential  V  .  Since  the  annular  velocity  vector  u  is 
assumed  constant,  the  rotation  axis  has  an  invariable  direction 
in  space,  which  we  may  take  as  the  z-axis,  so  that  (1-6)  holds. 
Then  the  centrifugal  force  wx  (wx  x)  may  be  expressed  as  the 
gradient  of  the  centrifugal  potential  (1-8)  : 

*  *  Y  n2(x2  + y2)  ,  (6-5) 

and  (6-2)  becomes 


p(u  +  2w  x  u)  .  p  grad  W  +  di  v  P  ,  (6-6) 

the  gravity  potentia'  W  being  defined  by 

W  *  V  +  e  (6-7) 

as  usual . 

The  equilibrium  state.  For  equilibrium  we  have  w  0 
since  there  is  no  relative  movement.  Denote  the  (in  general  variable) 
equilibrium  density  by  oq  «  and  the  correspondi ng  potential  by 
W0  .  The  stress  tensor  Pr,  is  assimed  to  correspond  to  hydrostatic 
equilibrium,  so  that 

£o  -  -Pol  (6-8) 


1. 1  iii,riiirtfiMgttaiiliirifiv\r . . 


by  (5-16);  po  is  also  variable  in  general. 

Then  (6-6)  reduces  to 

0  *  p0gradW0  +  div£0  ,  (6-9) 

or  by  (6-8) 

0  *  pogradklo  -  grad  po  ,  (6-10) 

which  is  again  the  well-known  condition  for  hydrostatic  equilib¬ 
rium  (2-27). 

Deviations  from  equilibrium.  Let  a  material  particle  in 


B 

x 

A  u 

Or  — 

*0 


FIGURE  6.1.  Displacement  from  equilibrium 

the  equilibrium  state  be  at  point  A  ,  having  the  position  vector 
x  ,  and  let  the  same  material  particle  after  deformation  be  at 
point  B  ,  having  the  position  vector  £  i  cf.  Fig*  6.1.  All  this 
is  quite  analogous  to  sec.  2;  cf.  Fig.  2.1. 
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Now  the  left-hand  side  of  (6-6)  is  small  together  with  u^  , 
and  so  is  the  right-hand  side 

R  *  pgrad  W  +  div  P  ,  (6-11) 

so  that  it  is  irrelevant,  to  the  first  order,  whether  we  take  R 
at  A  or  B  .  On  the  other  hand,  the  single  terms  pgrad W  i ’d 
div£  are  not  small,  so  that  both  must  be  taken  at  A  or  both 
at  B  ,  but  consistently.  We  take  them  at  point  B  .  Thus 

-  ’  -b  “  pB9rad  WB  +  div-B’  (6-12) 

Now 

*o  *  7  n2(*5  +  y2Q)  » 

$  *  |  n2(x2  +  y2 ) 

differ  only  by  negl igi bly  smal 1  quantities  since  the  centrifugal 
potential  itself  is  already  small.  Thus  we  may  put  &0  *  *  ,  so 

that  the  potential  disturbance  is  purely  gravitational: 

W  -  W0  =  (V  +  ♦)  -  (V0  +  ♦  )  -  V  -  V0 

-  Ve  +  Vi  ,  (6-13) 

consisting  of  the  tidal  potential  V  and  the  change  Vj  of 
the  body's  attractional  potential  caused  by  its  change  of  shape 
due  to  the  deformation;  cf.  (2-22). 
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We  again  denote  the  difference  WB  -  WA  by  n  : 

n  =  WB  -  WA  »  u_  •  gradW  *  _u  •  grad  W  o  = 

=  W0B  -  W0A  ,  (6-14) 

by  (2-33).  Thus  (6-13)  gives 

WB  *  WoB  +  (Ve  +  Vi  )B  , 
or 

WB  =  W0A  +  n  +  Ve  +  Vj  ,  (6-15) 

using  (6-14)  and  the  fact  that  Ve  and  Vi  are  so  small  that  it 
is  irrelevant  whether  we  take  them  at  A  or  B  . 

Let  us  now  find  a  similar  expression  for  pb  .  During  defor¬ 
mation,  a  mass  element  dM  may  change  its  shape  but  not  its  mass. 
Thus , 

dMB  »  dM0A  .  (6-16) 

On  the  other  hand,  the  corres pondi ng  volume  element  dv  changes 
by 

dvB  *  (1  +  ®)dv0A  (6-17) 

by  (5-8) ,  where 


Q  «  di v  u  . 


(6-18) 


Since  dM  *  pdv  ,  we  have  from  (6-16)  : 

pBdvB  B  poa^oa  » 
and  by  (6-17) 

pB ( 1  +  0 )  *  p0A  (6-19) 

or 

Pb  s  poa  "  PoQ  •  (6-20) 

Since  0  is  very  small,  p0  in  the  last  term  may  refer  to  A 
or  to  B  . 

Let  us  now  put 

Pi  =  PB  "  p  Ob  a  p  i  b  °1A  *  PA  “  POA  (6-21) 

for  similar  reasons  as  for  (6-15).  Then  (6-20)  gives 

Pi  a  pB  -  p0B  =  p0A  *  PoQ  ~  pob  •  (6-22) 

Now , 

PoB  *  poa  +  1  *  9radPo  (6-23) 

by  Taylor's  theorem.  The  combination  of  (6-22)  and  (6-23)  final 
ly  gives 

Pl  a  -p0o  -  o  *  gradp0  .  (6-24) 


We  also  note  that  the  change  of  the  gravitational 

potential  is  produced  by  the  density  anomaly  i  and  is, 
therefore,  related  to  it  by  Poisson's  equation  (Heiskanen 
and  Moritz,  1967,  p.5) 

AV!  *  -4irGpi  *  4irGdi  v(poju )  (6-30) 

by  (6-25).  Here  a  denotes  the  Laplacian  operator  and  G  the 
gravitational  constant. 

Let  us  finally  turn  to  the  stress  tensor  P  .  We  define  the 
incremental  stress  tensor  £i  by 

Pi  -  £b  -  Poa  •  (6-31) 

It  thus  represents  a  Lagrangian  increment  since  P  and  Pp  re" 
fer  to  the  same  material  particle  but  to  different  spatial  points. 
This  is  quite  analogous  to  (2-32), 

Now  we  can  substitute  (6-15),  (6-20),  and  (6-31)  into  (6-12): 

R  *  (pqa  “  po0)grad(WOA  +  n  +  Ve  +  Vi)  +  div(P0A  +  Pi ) 


=  O0A9radW0A  +  divPoA  " 


-p09gradW0  +  p0grad(n  +  V@  +  Vi)  +  divPi 

We  now  take  (6-9)  at  the  point  A  and  subtract  it  from  the  last 
equation.  There  remains 

R  =  -p0OgradW0  +  p0grad(n  +  Ve  +  +  div£i  .  (6-32) 


.  .  ►aM t.-  -  ■ 


jfctrfo  Wtfoaa: .  .mk.  - ? 
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By  the  definition  (6-11),  this  is  the  right-hand  side  of 
the  equation  of  motion  (6-6),  so  that  this  equation  becomes 

•«  * 

+  2u  x  jj)  »  -pQgradW  +  pgrad(n  +  V*  +  Vi  )  +  div£i  . 

(6-33) 

% 

Since  now  all  summands  are  small  of  first  order,  we  have  been 
able  to  replace  pq  and  Wo  by  p  and  W  ,  apart  from  second- 
order  quantities.  This  simplifies  the  notation;  still  it  is  usual 
ly  convenient  to  refer  p  and  W  in  (6-33)  to  the  equilibrium 
state . 

Eq.  (6-33)  represents  the  basic  equation  of  motion  for  our 
present  purposes. 

The  elastic  stress  tensor  Pi.  We  note  that  Pi  represents 

the  elastic  part  of  the  stress  tensor  £  .  In  fact,  £  consists 
of  a  large  hydrostatic  part  £o  which  provides  the  huge  pressure 
necessary  to  resist  the  gravitational  compression,  and  a  small 
residual  part  Pi  which  may  be  considered  to  be  related  to  the 
small  elastic  deviations  u^  by  Hooke's  law  (5-17)  or  (5-18). 

Thus , 

Pi  *  xel_  +  2yE  .  (6-34) 

Putting 

Pi  =  [Pij]  ,  (6-35) 


we  may  write 
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P 


U 


.*tj 


(6-36) 


in  view  of  (5-4)  and  (5-17). 

Various  forms  of  the  basic  equation.  The  right-hand  side 
R  of  the  basic  equation  of  motion  (6-33)  may  be  written  in  differ 
ent  forms.  In  this  way  it  is  seen  that  the  motion  equations  used 
by  different  authors  are  equivalent  in  spite  of  their  different 
forms . 

It  is  convenient  to  use  index  notation.  Then  (6-32)  becomes, 
on  again  omitting  the  subscript  zero  and  using  f 6 - 14 )  and  (6-26), 


3  u.  3  W  3  U.  3W  32W 

R  ■  -p  — —  -  +  p  - +  P  U.  -  + 

3xt  3X1  3xfc  3x.3xr 

3  3P.. 

♦  0 - (V,  +  V,)  ♦  — il  .  (6-37) 

3  X  .  3  X  . 

x  J 


This  form  has  been  qiven  by  Jeffreys  and  Vicente  (1957a,  p.  144) 
and  used  by  Wahr  (1981a).  Smith  (1974,  p.  494)  has 


R. 


Ifu. 

3 

_  9U/ 
p0  - 

♦  — 

Po 

3  X  ^ 

3Xi 

l  »Xkj 

3xk 

l 

3xi. 

3 

3*W 

+  p  -  (Ve 

+  Vi)  + 

PU  - 

I 

3  X 

i 

3Xi 

axk 

(6-38) 
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where  po  is  the  hydrostatic  pressure  in  (6-8).  By  differen¬ 
tiating  the  products,  this  is  readily  brought  into  the  form 


R 


i 


3p  3p0  3u.  apo  au  a 

- - *  + - £•  +  P  -  (Ve  +  Vi)  + 

ax.  ax,  ax,  ax,  ax.  ax. 

3  i  k  k  i  i 


a2w 


+  pu. 


ax,  ax, 

i  k 


Since,  by  (6-10), 


(6-39) 


9  Po  a W0 

-  a  p0  -  ,  (6-40) 

3xi  3Xi 

this  is  immediately  seen  to  be  equivalent  to  (6-37). 

The  concise  form  given  by  Dahlen  (1972,  p.  360): 


aw  a  a  ap.  . 

r  *  pi  -  +  -  (Pn)  +  p  -  (V  +  V.  )  +  — ii  (6-41) 

a  x .  ax.  ax.  ax, 

ii  i  d 

can  be  transformed  to  (6-37)  in  a  similar  way. 

Molodensky's  equation.  Let  us  return  to  the  form  (6-33)  : 

u  +  2w  *  u  =  grad(Ve  +  Vi  +  n)  -  SgradW  +  p" 1  di  v£i  .  (6-42) 

It  presupposes  that  the  z-axis  coincides  with  the  rotation  axis 

and  that  the  rotation  is  uniform.  In  other  terms,  it  holds  for  the 
nutation  frame  (sec.  1). 


The  corresponding  equation  of  motion  in  th i  body  frame 
is  found  as  in  sec.  1  :  we  must  add  (1-34)  to  tne  left-hand 
side,  or  subtract  it  from  the  right-hand  si de .  The  result  is 

x  u  =  grad  +  Vi  +  $  +  n  i  - 

"  2  £  H  i3  -  QgradW  +  i  divP!  ,  (6-43) 

where  w  *  (0,0, a)  as  throughout  the  present  section.  This 
equation  has  been  given  by  Molodensky  (1961,  eq .(11));  see  also 
(Jobert,  1964,  p.  66). 

An  incompressible  fluid  may  formally  be  considered  a  special 
case  of  an  elastic  medium  with  0=0  and 

?i  =  -PiJL  ,  (6-44) 

Pj  being  an  incremental  hydrostatic  pressure,  as  we  have  seen 
at  the  end  of  the  preceding  section.  For  this  particular  case, 
(6-43)  reduces  to  (2-1)  with  (2-38)  and  (2-39),  as  it  should. 


DIFFERENTIAL  EQUATIONS  FOR  THE  MANTLE 


Molodensky's  equation  (6-43)  holds  for  a  rotating  earth 
of  arbitrary  shape,  in  particular  for  an  ellipsoidal  earth.  It 
furthermore  holds  for  an  elastic  and  for  a  liquid  body,  a  liquid 
being  formally  considered  a  limiting  case  of  an  elastic  body  as 
we  have  seen  at  the  end  of  sec.  5. 

Eq .  (6-43)  also  holds  for  a  body  that  is  partly  elastic 
and  partly  liquid,  such  as  for  an  earth  consisting  of  an  elastic 
mantle,  a  liquid  outer  core,  and  a  solid  elastic  inner  core,  pro¬ 
vided  we  add  appropriate  boundary  conditions  at  the  boundary  sur¬ 
faces  between  core  and  mantle  and  between  outer  and  inner  core. 

Thus  (6-43)  may  be  considered  a  system  of  partial  differen¬ 
tial  equations,  to  be  solved  for  the  displacement  vector  jj  and 
for  other  quantities.  Such  a  uniform  treatment  of  (6-43),  together 
with  the  boundary  conditions  mentioned,  for  the  whole  earth  is 
particularly  suited  for  the  use  of  a  modern  computer;  it  has  been 
employed  by  Shen  and  Mansinha  (1976)  and,  most  recently  and  satis¬ 
factorily,  by  Wahr  (1979);  see  sec.  li. 

Molodensky  (1961),  however,  persued  a  largely  analytic  ap¬ 
proach,  as  did  Jeffreys  and  Vicente  (  1957a, b)  before  him.  Here  it 
proves  appropriate  to  treat  core  and  mantle  in  a  completely  separa¬ 
te  way.  The  core  is  treated  by  methods  described  in  Part  A  of  the 
present  report,  whereas  for  the  mantle  eqs .  (6-43)  are  used,  but 
in  a  greatly  simplified  form. 

Spherical  approximation.  In  fact,  the  consideration  of  ellip 
ticity  is  essential,  above  all,  for  the  core  as  we  have  seen  in 
sec.  2.  For  the  mantle,  a  spherical  approximation  is  sufficient  as 
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has  been  recognized  already  by  Jeffreys  (1949). 

Molodensky  neglects  the  flattening,  the  rotational  effects, 
and  even  the  free  elastic  vibrations  of  the  mantle.  He  thus  puts, 
formal ly , 


8  *  0  ,  *  =  0  ,  u  =  0  (7-1) 

(only  for  the  mantle).  These  neglections,  surprising  as  they  are 
at  first  sight,  can  be  justified  by  estimating  orders  of  magnitude 
and,  above  all,  by  good  agreement  with  the  more  precise  theory.  An 
essential  simplification  is  achieved  in  this  way. 

Thus  (6-43)  reduces  to 

0  =  pgrad(Ve  +  Vi  +  n )  -  pQgradW  +  divPi  .  (7-2) 

In  the  sequel  we  shall  more  or  less  follow  (Jobert,  1964). 

Spherical  coordinates.  For  a  spherical  earth  it  is  appro¬ 
priate  to  use  spherical  coordinates  r  (radius  vector),  e  (polar 
distance),  and  x  (longitude);  cf.  (2-41). 

In  order  to  express  (7-2)  in  terms  of  r,e,x,  we  first  note 
that  the  gradient  of  a  function  F  has,  in  these  coordinates,  the 
components 


gradF 


if. 

3r 

1  dF 
r  39 

1  dF 
rsine  3 X 


(7-3) 


.... - ,  . 


whereas  in  rectangular  coordinates  it  has  the  components 


gradF  = 


(7-4) 


More  difficult  is  the  expression  of  the  components  of  the  stress 
tensor  £1  in  spherical  coordinates.  From  standard  texts  (e.g., 
Love,  1927,  pp.  56,  91;  McConnell,  1931,  p.  311)  we  get  for  the 
stress-strain  equations  (6-36)  in  coordinates  r,6,x  : 


3U 

I©  +  2y  - - 

3  r 


X  0  +  2  u - - 

r  30 


fl  ^u  u  ] 

! - 1  +  JL\ 


=  X0  +  2y 


1  3ux  ue 

- -  +  -1  cote  + 

rsine  3X  r 


(7-5) 


U  f  1  3Ur  3Ua 

r  (  sine  3 x  3r 


"  ux  j  * 
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pex 


W  3U  1  3U  1 

-  — -  + - u.cote  ! 

r  [  36  sine  3X  j 


(7-5) 


Here  u^,  uQ  ,  u^  denote  the  spherical  components  of  the 
displacement  vector  u  ,  and  similarly  prr,  prQ  ,  etc.  denote 
the  components  of  the  stress  tensor  .  The  Lame  parameters, 
formerly  designated  by  X,  u,  are  not  denoted  by  T,  p,  in 
order  to  avoid  confusion  with  longitude  x 

These  rather  complicated-looking  equations  become  immediate¬ 
ly  clear  if  we  write  (6-36)  in  rectangular  coordinates  xyz  : 
putting  pxx  =  pu  ,  pxy  =  p 1 2  ,  etc.,  we  have 


3  U 

p  *  X0  +  2u  — -  , 


etc . 


Thus  there  correspond 


3  x ,  3y ,  3  z  and  3r,  r36,  rsin03X  , 


(7-6) 


(7-7) 


in  agreement  with  (7-3)  and  (7-4).  This  fully  explains  the  terms 
in  (7-5)  containing  partial  derivatives;  remaining  terms  containing 
ur  ,  uQ  ,  u^  directly  are  due  to  the  curvilinear  character  of 
the  spherical  coordinate  system. 

For  the  spherical  components  f  ,  f0 ,  fx  of  the  elastic 
force  vector 


vi’ u  n  r-  if  i  ■ 


f  =  di  v  £  i 
we  get 

3p 


(7-8) 


*P, 


„  rr  ,  '  r0 

rf  «  r  -  +  -  + 

r  3  r  30 


3pre  ,  5pee 

rU  * r  —  — 


3  r 
+  3p 


36 


r0 


3PrX  ,  3P8X 
rf  *  r -  +  -  + 


3  r 


1 

3prX 

s  i  ne 

3  X 

1 

3pex 

si  ne 

3  X 

1 

3pxx 

sine 

3  X 

+  2pcc- P98  -  Pu  +  precot9 

+  tpee  -  pxx)i:°te  + 


(7-9) 


*  3ptx  *  2pexc°t9  • 


Again  the  analogy  to  the  rectangular  case 


sp  3  p  3  p 

f  .  +  _**.  +  ._*£  ,  etc.  (7-10) 

3  X  3y  32 

helps  to  make  the  spherical  equations  transparent. 

The  volume  dilatation  (5-11)  takes  the  form 

3  u  1  3u  „  1  3u .  2  cote 

e-  —  + - -  + - -  +  -  ur  +  -  u0  .  (7-11] 

a  r  rae  rsineax  r  r 
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Second-degree  harmonic  displacements.  In  order  to  simplify 
and  solve  the  system  of  partial  differential  equations  (7-2),  it 
is  customary  to  take  for  the  components  of  the  displacement  vector 
u  the  following  special  form: 


ur  a  H  (  r ) S  ( Q  ,A  )  =  HS  , 

“6  *  L<r!  *  Lr“se  •  <7-17> 

ui  *  rsi ne  TT  *  Mrsin«)  1  Sx  • 

where  SQ  and  denote  partial  derivatives  (in  contrast  to 

u0  ,  ux  which  simply  are  components). 

Here  H(r)  and  L(r)  are  functions  of  the  radius  vector  r  , 
whereas  S(0,a)  denotes  any  Laplace  surface  harmonic  of  second 
degree.  For  the  present  purpose  (tidal  effects  on  the  direction 
of  the  rotation  axis)  it  is  appropriate  to  consider  S(e,A)  a 
spherical  harmonic  of  degree  n  »  2  and  order  m  =  1  ,  that  is, 
of  the  form 


S ( e , x )  *  P11(cos0)(acosA  +  bsinA) 
*  aRj i ( e  ,  a  )  +  bSi i ( 0 ,  a  )  , 


(7-18) 


where  a  and  b  are  arbitrary  constants.  For  notations  cf. 
(Heiskanen  and  Moritz,  1967,  p.  29),  and  for  a  justification  of 
the  restriction  to  n  =  2  cf.  (Moritz,  1980a,  sec.  55);  note  also 
that  (7-18)  is  in  agreement  with  the  form  (2-40)  of  the  present 
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report  for  the  tidal  potential. 

Being  a  Laplace  harmonic  for  n  *  2  ,  the  ^unction  S(e,x) 
satisfies  the  partial  differential  equation 

see  +  sTRTT  saa  *  'ot6Se  -  -  6S  <7-19> 

4 

(Heiskanen  and  Moritz,  1967,  p.  20,  eq.  (1-44)). 

We  also  note  that  Molodensky  (19G1)  uses  the  notation  T(r) 
instead  of  L(r)  ;  we  have  restricted  the  letter  T  to  toroidal 
oscillations,  so  that  our  notation  L(r)  is  consistent  with 
sec.  3. 

With  (7-17)  let  us  now  compute  the  volume  dilation.  We  find 

au 

—  *  H  '  (r)S(e  ,X  )  *  H'S  ;  (7-20) 

ar 

primes  will  always  denote  d/dr  ,  so  that  H'  =  dH/dr  .  Similarly 


a  u  au, 

— -  -  Lr‘iSao  ,  — - 

ae  9  ax 


L(rsine  J"1  S 


xx 


(7-21) 


Thus  (7-11)  becomes 

O  .  H'S  *  2r-‘HS  *  r-*L(S„  ♦  ^  SJA  *  cotes,) 
»  (H‘  +  2r-1H  -  6r“2L )S 

by  (7-19).  This  may  be  written  in  the  form 
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0  *  XS  (7-22) 

wi  th 

X  *  H'  +  2r*lH  -  6r"2  L  .  (7-23) 

In  the  same  straightforward  way  we  get  from  (7-5)  the  comp 
nents  of  the  stress  tensor  : 

Prr  *  NS  , 

Pee  =  +  ^ur  2 LSee  • 

Pxx  *  QS  +  2ur"2 L(si n’2eSxx  +  cot6S0)  , 

(7-24) 

pr9  *  --3NS6  ■ 

PrX  a  r*3Msin’1eSx  , 
pQX  *  2wr"2Lsin~16 (S0X  -  cot0Sx) 
with  the  abbreviations 

M  *  u  ( r2  H  +  r2 L '  -  2rL)  ,  (7-25) 

N  «  XX  +  2uH 1  ,  (7-26) 


Q  «  “X  +  2u  r" 1 H 


(7-27) 
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Note  that  N  is  proportional  to  the  normal  tension  prr  , 

'and  that  M  is  proportional  to  the  tangential  tensions  prQ 
and  prX  acting  on  a  surface  r  ■  const. 

In  the  same  way  we  get  for  the  components  of  the  elastic 
force  (7-9)  : 

fr  ■  (N1  -  6r~4M  +  4ur“ 1 H 1  -  4pr“2H  +  12wr’3L)S  , 

fQ  *  ( r“ 3M 1  +  r-iQ  -  lOur-3L)S0  ,  (7-28) 

fx  ■  (r ** 3 M 1  +  r_1Q  -  lOy r“3L )si n" 1 0SX  . 

Once  more  we  remark  that  the  subscripts  r,e,x  in  prr,  prQ,  ... 
....  fx  denote  components,  whereas  in  S0,  S0X,  ...  they  denote 
partial  derivatives. 

The  quantity  n  is  immediately  obtained  from  (7-15)  and  (7-17) 

n  *  HSW  *  .  (7-29) 

We  finally  try  the  representation 

Ve  +  V !  ■  R( r )$( 0 , A )  =  RS  .  (7-30) 

A  first  justification  is  provided  by  the  fact  that  the  external 
tidal  potential  V  is  indeed  a  spherical  harmonic  with  n  *  2  , 
m  =  1  ,  and  that  the  response  V5  to  the  disturbance  can 

be  expected  to  be  proportional  to  S(e,x)  as  well.  The  final 
justification  lies  in  the  fact  that  the  substitution  (7-30)  works, 
as  will  be  seen . 
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Now  we  can  write  (7-2)  In  spherical  coordinates: 
fr  +  p  Tr  +  V1  4  n  )  '  P0W‘  *  °  ’ 

fQ  +  p  r’1  Tfr  (ve  +  +  n)  *  0  ,  (7-31) 

fx  +  p(rsine)'1^-  (Vfi  +  Vi  +  n)»  0  ; 

here  we  have  used  (7-3),  (7-8),  and  (7-13).  We  substitute  (7-28), 
(7-29),  and  (7-30),  differentiating  according  to 

■ft  (Ve  +  ♦  l)  *  S  £  (R  +  HU')  • 

•  S( R 1  +  H'W'  +  HW “ )  . 

^  (Ve  +V1+n)*(R+  HW')Se  . 

The  first  two  equations  of  (7-31)  then  give 

N'  -  6r-'*M  +  4yr“ 1 H '  -  4yr"2H  +  12ur"3L  + 

+•  p(R'  +  H'W'  +  HWM  -  XW1)  -  0  ,  (7-3?) 

r" 3M '  +  r’lQ  -  lOyr"3!.  ♦  r’lp(R  +  HW')  »  0  ;  (7-33) 

the  third  equation  of  (7-31)  provides  a  check  for  (7-33). 


’^‘rn'prTTT^  c*,ri  v? ' 
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This  is  the  result  of  the  el  astici  ty equati ons  (7-2).  They 
are  not  yet  sufficient  to  fully  determine  the  problem.  We  still 
need  Poisson’s  equation  (6-30)  , 

a •  -4if 6p i  ■ 

*  4itG(u  *  gradp  +  p©) 

by  (6-24),  writing  p  for  pq  as  usual.  For  spherical  symmetry, 


o  ■  p  ( r)  , 


this  becomes 


(7-34) 


AVi  »  4irG(p'ur  +  p<3)  *  4nG(p'H  +  pX)S 


(7-35) 


Since  the  tidal  potential  Ve  is  harmonic  (outsideof  sun  and  moon) 


we  have 


aV  *  0  f 

e 


(7-36) 


so  that 


4{Vft  +  V:)  *  A(RS)  =  A V i 


(7-37) 


as  given  by  (7-35).  From  (7-16),  (7-19),  and  (7-30)  we  obtain 


A ( RS )  =  (R"  +  2r*lR'  -  5r*2R)S  . 


(7-38) 


The  comparison  of  (7-35)  and  (7-38),  using  (7-37),  finally  gives 


R"  +  2r_1 R '  -  6r“2  R  *  4tt  G  (p  *  H  +  pX)  .  (7-39) 

We  now  introduce  a  new  auxiliary  variable  P(r)  by 

P  *  r2  ( R '  -  4tr  Gp  H )  (7-40) 

(Molodensky  denotes  it  by  L  ).  We  solve  (7-40)  for  R'  ,  differ¬ 
entiate  to  obtain  R"  and  substitute  into  (7-39).  The  result  is 

P'  =  6R  +  4ttGp  r2  ( X  -  H'  -  2r-1H)  .  (7-41) 

The  system  of  the  two  first-order  differential  equations 
(7-40)  and(7-41)  is  equivalent  to  the  second-order  differential 
equation  (7-39);  the  advantage  is  that  now  p'  is  eliminated, 
so  that  p  need  not  be  differentiated.  This  is  important  since 
the  differentiation  of  an  empirically  (and  not  too  well)  determi¬ 
ned  quantity  such  as  p  is  problematical. 

Equations  (7-25),  (7-26),  (7-32),  (7-33),  (7-40),  and  (7-41) 
will  be  our  basic  equations.  After  elimination  of  X  and  Q  by 
(7-23)  and  (7-27)  they  become 

M  =  y  (  r2 1 '  -  2rL  +  r2H)  , 

N  *  (X  +  2w )H '  +  A ( 2 r" 1 H  -  6r'2l)  , 


G 
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(7-42) 

P'  *  6R  -  24*GpL  , 

N 1  -  6r"4M  +  4ur" 1 H '  -  4yr‘2H  +  12yr-3L  i- 
+  p ( R '  +  HW"  -  2r“ 1 HW 1  -  6r‘2LW')  =  0  , 

r- 3  M '  +  Xr" 1 H  *  +  2{x  +  p)r"2H  -  (6x  +  10u)r"3L  + 

+  pr*1  (R  +  HW  )  =0  . 

This  is  a  system  of  6  ordinary  differential  equations  of  the  first 
order  for  the  6  unknown  functions  H(r),  L(r),  M (r),  N(r),  P(r), 
R(r)  .  The  density  p(r)  ,  the  lamfe  parameters  x(r)  ,  p(r)  (we 
have  again  written  X  instead  of  X  )  ,  and  the  potential  W(r) 
are  considered  known  from  a  standard  earth  model. 

A  basic  advantage  of  the  system  in  its  present  form,  which 
goes  back  to  (Molodensky,  1953),  is  that  it  does  not  contain  deriv¬ 
atives  of  the  empirical  quantities  x,y,p  . 

The  equation  system  (7-42)  will  be  further  considered  in  the 
next  section. 


I 
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8.  SOLUTION  FOR  THE  MANTLE 

The  basic  equation  system  (7-42)  may  be  solved  algebraically 
for  the  derivatives  H',  L',  R'  .  The  result  has  the  form 

H1  =  a  1 1  H  +  ai  2  L  +  a  1 3  R  +  a  ^  4  M  +  a  1 5  N  +  a  1  g  P  , 

L'  =  32  1  H  +  322^-  +  323^  +  324  M  +  32sN  +  326  ^  » 

R'  *  a  3 1 H  +  332  L  +  833R  +  a34M  +  a3sN  +  a3gP  , 

(8-1) 

M '  =  84 1 H  +  a  4  2  L  +  3438  +  a44M  +  345  N  +  a  4  g  P  , 

N1  =  as  1 H  +  as2 L  +  as3R  +  as4M  +  assN  +  asgP  , 

P‘  =  ag  i  H  +  a&2  L  +  *6  3  R  +  ag4M  +  36sN  +  ae6P  • 

The  coefficients  aiR  are  functions  of  r  .  For  instance,  the 
comparison  between  the  last  equation  of  (8-1)  and  the  fourth 
equation  of  (7-42)  shows  that 

36  1  =  0  *  a62  =  ~2 4 -it Gp  ,  a6  3  =  6  , 

36  4  =  a65  =  3  g  6  =  0  • 

If  x.w.p,  W  are  given  as  functions  of  r  ,  then  the  coefficients 
aik(r)  can  be  computed  and  the  system  of  differential  equations 
(8-1)  integrated  by  a  numerical  integration  technique,  provided 
suitable  boundary  values 

Ho »  Lc »  Ro >  Mq »  N0 ,  Po »  (8-2) 

are  given.  We  take  them  to  be  the  values  of  the  corresponding 
quantities  at  the  earth's  surface,  considered  a  sphere  of  radius 
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a  »  6371  km  : 

H0  -  H(a ) ,  L0  «  L(a),  . ,  P0  =  P(a)  .  (8-3) 

A  slightly  different  approach  is  the  following.  We  solve 
the  system  (8-1)  six  times,  first  with  boundary  values 

H0  *  1 ,  Lo  a  Ro  a  *  ’  •  =  Po  =  0  »  (8-4) 

denoting  the  solution  by  Hi(r),  Li(r),  Pi(r)  .  Then  we 

solve  with  boundary  values 

H0  -  0,  L0  •  1,  R0  -  M0  »  N0  -  P0  -  0  ,  (8-5) 

obtaining  solutions  H2(r),  L2(r),  . .  P2(r)  .  We  continue 

with  Rq  *  1  ,  all  other  boundary  values  being  zero,  etc.  Finally 
we  take 


Ho  s  Lq  *  ■  '  •  s  Nq  s  0  ,  Po  =  1  ,  (8-6) 

denoting  the  solution  by  H0(r),  L6(r) . P  6  ( »" )  • 

Then  the  other  solution  of  (8-1)  for  general  boundary  values 
(8-2)  is  given  by 

H(r)  *  H0Hi(r)  +  L0H2(r)  +  •  •  •  +  P0H6(r)  , 

L(r)  =  H  o  L  i  (  r )  +  LoL2(r)  +  *  *  *  +  Po  ^-6  (  r )  , 

(8-7) 


P ( r)  =  H o P i ( r )  +  L0P2(r)  +  •  •  •  +  P0P6(r)  . 
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This  linear  dependence  of  the  solution  on  the  boundary  values  fol¬ 
lows  from  the  linearity  of  the  differential  equations  (8-1):  a 
linear  combination  of  solutions  is  also  a  solution.  The  boundary 
conditions  are  also  satisfied  since,  for  r  »  a,  (8-7)  gives 

H(a )  *  Ho  •  1  +  Lo  •  0  +  •  •  •  +  Po  •  0  *  Ho  ,  etc. 

These  two  approaches,  the  direct  approach  using  the  boundary 
values  (8-2),  and  the  indirect  approach  with  (8-4),  (8-5),  .... 
correspond  to  the  solution  of  a  system  of  algebraic  linear  equa¬ 
tions,  first  by  a  direct  method  such  as  Cholesky's,  and  second 
by  means  of  the  inversion  of  the  matrix. 

The  second  approach  is  particularly  transparent  and  simple 
in  the  present  case,  in  view  of  the  special  form  of  the  boundary 
values  to  be  considered  now. 

The  boundary  values.  First  we  note  that  two  boundary  values 
are  even  zero.  This  follows  from  the  condition  that  the  earth's 
surface  is  free  from  stress.  The  stress  vector  £  acting  on  the 
sphere  r  *  a  is  given  by  (5-13),  with 

n  -  (1,0,0) 

for  the  radial  unit  vector  in  spherical  coordinates  (only  the  first, 
radial,  component  of  £  is  different  from  zero).  We  get 
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which  must  be  zero  at  the  earth's  surface  r  *  a  .  Thus  all 
tensions  prr  (normal)  and  pr0  and  prX  (tangential)  vanish 
for  r  ■  a  .By  (7-24)  this  leads  to  the  conditions 


M0  -  M ( a )  -  0  ,  (8-9a) 

N0  -  N(a )  *  0  .  ( 8-9b ) 

We  further  put 

H0  -  H(a)  *  h  ,  (8-10) 

L0  -  L ( a )  »  1  .  (8-11) 


For  the  external  potential  Ve  we  have  fojnd  the  expres¬ 
sion  (2-43)  : 
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»c  r2Si 


(8-12) 


where 


Si  «  Si(9,x,t)  »  si ne cosecos (o t  -  x )  .  (8-13) 

For  r  *  a  this  gives 

Ve(a)  -  <a2S1  ,  (8-14) 

so  that,  by  (8-12) 

V4(r)  ■  £  V#( a)  .  (8-15) 


1 

3 


I 
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It  is  clear  that  (8-14)  represents  a  spherical  surface 
harmonic  (n  ■  2,  m  «  1),  and  we  may  therefore  define  the  function 
S(e  ,x )  in  (7-17)  to  be  identical  to  IT  (a)  : 

S(0,x)  -  V#(a)  -  <  a2  Si  .  (8-16) 

Then  (7-30)  and  (8-15)  give 

a 

Vi(r)  -  R(r)S  -  Ve(r)  -  (R(r)  -  £)  S  (8-17) 

and 

Vi(a)  -  (  R(a )  -  1]S  «  [R(a)  -  1)V  (a)  .  (8-18) 

Denote  by 

k  *  R(a)  -  1  (8-19) 

the  factor  of  proporti onal i ty  between  the  indirect  effect  V i 
and  the  tidal  potential  at  r  *  a  .  It  is  nothing  else  than 

the  well-known  Love  number  for  the  potential. 

Thus , 


R0  *  R ( a )  »  1  +  k  (8-20) 

expresses  the  boundary  value  of  R  in  terms  of  the  Love  number  k. 
Now  we  see  that  also  h  and  1  in  (8-10)  and  (8-11)  can  be  identi¬ 
fied  with  the  usual  Love  numbers.  In  fact,  the  definition  of  the 
Love  numbers  (cf.  Moritz,  1980a,  p.  480)  reads 
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9 


(8-21) 


for  r  «  a  ,  where  g  is  mean  gravity  at  the  earth's  surface. 
Putting  r  »  a  in  (7-17)  and  comparing  to  (8-21)  we  find  with 


H ( a )  -  g“l  h  , 
L(a)  ■  ag_1l  . 


(8-22) 


If  we  chose  units  such  that 

a  *  1  ,  g  *  1  ,  (8-23) 

we  get  (8-10)  and  (8-11). 

The  last  boundary  value,  P0  ,  can  be  expressed  in  terms 
of  the  Love  number  k  by 

P0  =  P(a)  *  2  -  3k  .  (8-24) 

The  derivation  of  this  simple  relation  is  somewhat  lengthy  and 
will  be  deferred  to  the  end  of  the  present  section. 

Equations  (8-9)  through  (8-11),  as  well  as  (8-20)  and  (8-24), 
give  the  boundary  values.  Thus  the  solution  (8-7)  may  be  written 


n  'i  f ugaHSSTiii  liii  ) "  i:i 
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H(r) 

Mr) 


hHi  (r)  +  lH2(r) 
hLi (r)  +  ll2  (r) 


♦  (l  +  k)H3(r)  ♦  (2“  3k)H6(r) 
+  (l  +  k)L3(r)  ♦  (2-  3k)l6(r) 


I 


I 


P(r) 


(8-25) 


hPMr)  ♦  lPa(r)  +  (1  +  k)P3(r)  + 


(2  -  3k)P6(r)  . 


,The  boundary  values  on  which  the  solution  depends,  are  only 
the  Love  numbers  h,  k,  1  .  These  parameters  h,  k,  1  will  be 
determined  in  the  next  section  by  means  of  the  interaction  between 
core  and  mantle. 

The  functions  H4(r),  Hs(r),  L4(r),  Ls(r),  ....  P4(r),  Ps(r) 
do  not  enter  in  (8-25)  and  are  not  needed. 

Derivation  of  (8-24).  Consider  the  surface  of  the  earth 
before  and  after  deformation  (Fig.  8.1).  Consider  further  the 


after  deformation. 


potential  Vg  .  This  potential  may  be  regarded  as  being  produced 
by  the  small  density  distributed  inside  the  earth's  surface, 
combined  with  a  surface  layer  of  density 


-Vfo-  .1  •  .1  ni'-iii- 
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<  *  pur  .  (8-26) 

In  fact,  the  space  between  the  surfaces  z  and  s',  which  was 
empty  before  the  deformation,  is  filled  with  density  p  after 
the  deformation;  and  the  attraction  between  this  thin  layer 
between  i  and  z*  may  be  very  well  approximated  by  that  of  a 
material  surface  of  density  (8-26). 

Now  we  know  that  the  normal  derivative  3/3n  of  the  poten¬ 
tial  V  of  a  surface  layer  undergoes  a  discontinuity  at  the  sur¬ 
face 


f3V] 

3  Vl 

l  Tnj 

^  e  ^ 

[tttJ 

-4nGie  , 


(8-27) 


cf.  (Heiskanen  and  Moritz,  1967,  p.  6);  the  subscripts  e  and 
i  denote  exterior  and  interior  derivatives,  respectively.  In  our 
case,  the  normal  derivative  is  the  radial  derivative  3/3 r  ,  V 
is  Vj  ,  and  <  is  given  by  (8-26).  Thus, 


( 3  V  ^ 
l3r 


fav.il 

[3r 


'-4irGp  U. 


By  (8-17)  we  have 


(8-28) 


V i ( r )  -  [ R( r )  -  £  }  S  for  r  <  a  ,  (8-29) 

that  is,  inside  the  earth.  Differentiating  and  putting  r  =  a 
we  get 
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J|ij  « [R '  (a )  -  2a*1  ]S  .  (8-30) 

Outside  the  earth,  the  potential  V;  is  harmonic  as  any 
gravitational  potential  is.  Thus  it  is  given  by  the  expression 

Vi  (r)  *  [  R(a)  -  1  }*£  S(e  ,x )  for  r  >  a  .  (8-31) 

In  fact,  S(e,x)  is  a  spherical  surface  harmonic  of  degree 
n  *  2  ,  so  that  the  dependence  on  r  must  necessarily  have 
the  form 


for  a  harmonic  function  everywhere  regular  outside  a  sphere; 
cf.  (Heiskanen  and  Moritz,  1967),  pp.  20  and  34.  Furthermore, 
for  r  ■  a  it  must  reduce  to  (8-18).  These  two  conditions  unique¬ 
ly  determi ne  (8-31 ) . 

Differentiating  (8-31)  and  putting  r  »  a  we  get 

(|^j  -  -3a"1  [  R( a )  -  1 } S  .  (8-32) 

Now  we  substitute  (8-30)  and  (8-32),  together  with 

(u)  =  H(a)S  , 

r-a 


(8-33) 
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into  (8-28)  : 


-3a“l R(a )  +  3a"1  -  R'(a)  +  2a*1  »  -4»GpH(d) 


or 


R '  (a )  -  4nGpH(a)»  a-J[5  -  3R ( a ) ]  . 

By  (7-40),  the  left-hand  side  Is  a*2P(a)  ,  so  that 

P ( a )  -  a(  S  -  3R(a)  ]  . 

Substituting  (8-20)  we  find 
P(a)  -«  a ( 2  -  3k)  , 

and  with  (8-23)  we  get  (8-24),  which  was  to  be  derived. 
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9.  CORE  -  MANTLE  BOUNDARY  CONDITIONS 


As  in  the  last  section,  we  put 

a  •  1  (9-1) 

for  the  radius  of  the  earth's  surface,  approximated  by  a  sphere. 

To  the  same  approximation,  the  core-mantle  boundary  will  be  a 
sphere  of  radius 

b  «  0.55  .  (9-2) 

We  shall  now  derive  various  conditions  to  be  satisfied  on  the 
surface  r  «  b  . 

Parameter  M  :  no  tangential  tension.  The  comporents  of 
the  stress  vector  (8-8)  must  be  continuous  across  the  core-mantle 
boundary  r  «  b  .  Since  the  core  is  an  incompressible  fluid,  there 
is  no  tangential  tension  there.  By  continuity,  the  tangential  ten¬ 
sions  pr0  and  prA  must  be  zero  also  on  r  »  b  ,  which  by 
(7-24)  leads  to  the  condition 

M(b)  »  0  .  (9-3) 

By  (8-25)  this  may  be  written 

hM^b)  +  1 M2  ( b )  +  ( 1  +  k)M3  (b)  +  (2-3k)M*(b)  »  0  .  (9-4) 

Now  the  Mi(b)  are  constants  which  are  known  if  the  functions 
Mi(r)  have  been  computed  in  the  way  described  in  sec.  8.  There¬ 
fore,  (9-4)  gives  a  condition  to  be  satisfied  for  the  Love  numbers 
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h,  k,  1  .  Being  independent  of  the  core  structure,  it  is  the 
simplest  boundary  condition.  For  example,  Molodensky  (1961) 
gi  ves 


0.7053  k  -  0.4793  h  -  3.7155  1  +  0.3833  *  0  (9-5) 

for  his  Model  2  (cf.  sec.  10);  numerical  values  will,  of  course, 
be  slightly  different. 

Parameter  H  :  continuity  of  normal  displacement.  The 
normal  displacement  q  of  the  core-mantle  boundary  is  given 
by  (7-17)  : 

q  -  ur(b)  *  H(b)S(0,x)  .  (9-6) 

By  (7-15)  we  have  on  r  *  b 

n  *  qW'(b)  =  H(b)W'(b)S(9,x)  .  (9-7) 

Eq.  (8-16)  gives  with  (9-1) 

S(e,x)  **  k  Si  (9-8) 


where 


Si  =  si n9cos8 cos (o t  -  A )  . 


(9-9) 


In  the  liquid  core,  n  has  the  form  (2-47),  which  on  the  boundary 
r  ■  b  becomes 


[oj 


The  comparison  between  (9-7)  and  (9-10),  taking  account  of 
(9-8),  shows  that 

C  =  -<b-2H(b)W'(b)  .  (9-11) 

The  boundary  condition  expresses  the  parameter  C  in  (2-47)  in 
terms  of  H(b)  and,  by  (8-25)  for  r  =  b  ,  finally  in  terms  of 
the  Love  numbers  h,  k,  1  . 

Before  going  to  the  next  condition,  we  point  out  that  the 
spherical  approximation,  used  from  sec.  7  on,  is  legitimate  for 
the  mantle  including  its  lower  boundary  r  =  b  .  In  the  core, 
however,  the  ellipticity  plays  an  essential  role  and  must  be 
taken  into  account.  This  requires  some  care. 

Parameter  R  :  continuity  of  potential.  In  the  core  we 
have  by  (2-39)  with  (2-19)  and  (2-46)  : 

V  =  Ve  +  V!  +  (1  +  §)* 

a  (A  -  B)(xzco$ot  +  yzsi  not) 

■  (A  -  B)r2Si  .  ( 9 - 1 2  a  ) 

In  the  mantle  (7-30)  holds: 

Ve  +  V!  »  RS  .  ( 9- 12b  ) 

Since  the  potential  is  continuous  across  the  core  boundary, 
both  expressions  must  be  equal  for  r  «  b  : 


(A  -  B )b2Sl  =  R ( b ) < S x  . 


(9-13) 
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Here  we  have  used  (9-8)  and  neglected  the  term 


(1 


fl  +  o 

n 


(9-14) 


This  is  permissible  since  o  =  -n  for  the  main  tidal  frequencies 
and  since  the  quantity  4>  is  very  small.  Thus  (9-13)  gives 


A  -  B  =  <b'2R(b)  .  (9-15) 

This  equation  relates  the  core  parameters  A  and  B  to  R(b) 
and  hence,  by  (8-25)  to  h,  k,  1  . 

Parameter  N  :  equality  of  normal  stress.  Above  we  have 
pointed  out  the  continuity  of  the  stress  vector  across  the  boundary 
r  =  b  and  used  the  fact  to  derive  the  condition  (9-3).  Now  we 
take  the  normal  component  prr  .  In  the  core  this  corresponds  to 
■>Pi  where  pi  is  the  anomalous  pressure  (2-32).  Hence  prr  for 
the  mantle  and  -pi  for  the  core  must  be  equal  on  the  boundary 
r  =  b  ,  which  leads  to  the  condition 

Prr  +  Pi  s  0  .  (9-16) 

Now  prr  is  given  by  (7-24)  : 

Prr  a  NS  *  Ntc S i  (9-17) 

and  Pi  may  be  found  from  (2-43)  : 


Pi  *  p (A  -  C)r2Si 


(9-18) 
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in 


We  put  r  »  b  and  substitute  in  (9-16),  obtaining 

<N(b)  +  p b2  ( A  -  C)  «  0  .  (9-19) 

Here  we  may  eliminate  C  by  (9-11)  and  express  A  in  terms 
of  B  by  (9-15).  The  result  is 

pcolreN ( b )  +  R(b)  +  H(b)W'(b)  +  k “l  b2B  -  0  ,  (9-20) 

writing  p  =  pcore  for  the  (constant)  density  of  the  core. 

Parameter  P  :  discontinuity  of  the  potential  gradient. 

On  the  earth's  surface  r  »  a  -  1  we  have  the  condition  (8-28) 


-4rrGpur  , 


(9-21) 


expressing  the  discontinuity  of  aVi/ar  across  this  surface. 
This  is  due  to  the  density  jumping  from  p  inside  the  surface 
to  the  value  zero  outside,  so  that  p  in  (9-21)  represents 
this  density  jump. 

Across  the  core-mantle  boundary  r  ■  b  ,  the  density  jumps 

from  pcore  be1ow  t0  p mantle  above*  Thus  (9-21)  takes  for  the 
surface  r  «  b  the  form 


Sr  J 


mantle 


fay, 

1  3r  j 


'core 


■4wG( p  -p  ,)u  , 

'  core  mantle7  r 


(9-22) 


replacing  p  by  the  density  jump  o  -  o..  ;  the  notation 

cere  mantle 
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f®v,l 

m 

\p±] 

l  J 

e 

l  3r 

mantle 


is  clear  since  the  external  derivative  is  for  .id  by  approaching 
the  surface  r  «  b  from  above,  that  is,  through  the  mantle, 
and  similarly  for  (aV/3r)i  . 

We  may  write  (9-22)  in  the  form 


mantle 


4  it  Go  ,  U 

mantle  r 


which  shows  that  the  quantity 


(9-23) 


H 1  -  4itGpu  <9-24> 

3  r  r 

is  continuous  across  the  boundary  r  *  b  . 

The  external  tidal  potential  Ve  is  harmonic  inside  the 
earth  and  therefore  continuous,  together  with  its  derivative 
3V  / 3 r  ,  across  the  surface  r  *  b  .  The  same  holds  for  the 
incremental  centrifugal  potential  ♦  defined  by  (1-31).  Putting 

V  «  Ve  +  Vi  +  (1  ♦  §)* 

as  in  (9-12a)  we  therefore  conclude  that  the  quantity 


P-  -  4irGp  u 
3r  r 


(9-25) 


is  also  continuous  across  the  boundary  r  =  b  . 

In  the  mantle  we  may  neglect  <t>  and  replace  V  by  (9-lt.b), 
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so  that 


3  V 


3  r 


«  R'S  , 


3  V 


~  -  4*Gpur  =  (R*  -  4n Gp H ) S  «  Pr~2S  ; 


(9-26) 


we  have  used  (7-17)  and  (7-40). 

In  the  core,  V  is  given  by(9-12a),  so  that 


3  V 


jp  3  (A-  B)2rSi  , 


-4hGpu_  «  -3gr’ 1  u^  , 

r  r 


(9-27) 

(9-28) 


introducing  gravity  g  by  (2-61).  Thus 


|p  -  4ir Gp u^  *  (A-B)2r$i  -  3gr*lur  . 


(9-29) 


In  view  of  the  continuity  of  (9-25),  the  expressions  (9-26) 
and  (9-29)  must  be  equal  at  the  boundary  r  *  b  : 


P(b)b’2S  *  ( A  -  B ) 2bS i  +  3W  '  (b ) t>* 1 H( b )S 


Here  we  have  put  g(b)  •  - W ' ( b )  ,  so  that,  by  (2-61) 


(9-30) 


W'(b)  »  -  bp 


core 


(9-31) 


Finally  we  take  into  account  (9-15)  and  (9-8),  obtaining 


P(b)  «  2bR(b)  +  3bW‘(b)H(b) 


(9-32) 


iiv:  ^  w.tm 


vniilU. '  i  tk.lti,,!:-,1 , . 
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Equations  for  Love  numbers  and  fluid  parameters.  The 
basic  equations  are  the  boundary  conditions  (9-3),  (9-11),  (9-15), 
(9-20),  and  (9-32).  We  collect  them  now,  replacing  the  fluid 
cn;1s».«nts  A,B,C  by  the  equivalent  constants  o,b,y  related  to 
them  by  (2-20),  (2-64),  and  (2-66);  this  is  preferable  since 
the  symbols  A  and  B  will  be  used  for  the  principal  moments 
of  Inertia  and  since  a,8,y  are  dimensionless.  Thus  we  get, 
also  using  (9-31), 

a 

M  ( b )  ■  0  , 

e2b<"1Y“H(b)  *  0  , 

e2b-1W'  (b)*-1*-  0(o  +  2n)*-iB  -  b‘2R(b)  =  0  ,  (9-33) 

°^lr#N(b)  +  R(b)  +  H(b)W'(b)  +  b2o(o+ 2p)k' lB  »  0  , 

2bR(b)  -  P(b )  +  3bW' (b)H(b)  -  0  . 

We  combine  these  formulas  with  an  equation  following  from 
(2-58)  and  (2-65)  : 

7  (r-  a)e*  ♦  s  -  fe  =  0  .  (9-34) 

The  quantities  H(b),  M(b),  P(b),  R(b)  are  linear  combinations 
of  h,  k,  1  with  known  coefficients  H  i  ( b ) ,  H2(b),  etc.,  cf. 
(8-25);  the  constant  W'(b)  ,  given  by  (9-31),  is  likewise  known. 
Therefore,  (9-33)  and  (9-34)  represent  a  system  of  6  equations 
for  the  7  unknowns  h,  1,  o,  B,  y»  e. 


_N_  , ,  _  .  . _  . - _ ^ 
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Therefore,  If  we  could  find  an  independent  seventh  relation, 
then  we  could  compute  these  7  unknowns,  getting  a  unique  solu¬ 
tion  for  any  desired  value  of  the  frequency  o  .  Such  a  relation 
is  provided  by  the  Euler-Li ouvi 1 le  condition  to  be  discussed  in 
the  next  section. 


t 
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10.  THE  EULER  -  LIOUVILLE  CONDITION 


A  rigid  or  nonrlgld  body  rotating  with  an  angular  velocity 
vector  w  satisfies  the  equation 


+  «  *  H  -  L  ,  (1C-1) 

where  H  is  the  angular  momentum,  and  represents  the  torque 
(Moritz,  1980b,  p.  13);  this  is  Euler's  equation.  A  simple  con¬ 
sequence  is  Liouville's  equatior 


(Cu  +  h)  +u)  *  (C  w  +  h  )  *  JL 


(10-2) 


where  C  is  the  inertia  tensor  (ibid.,  p.  14),  and  hi  is  the 
relative  angular  momentum. 

We  put 

m  *  w_o  *  4  ,  C.  ■  C_q  +  £  ,  (10-3) 


0 

mi 

(jlQ  * 

1  0 

,  6w  ■ 

m2 

a  1 

! 

m3 

L 

L  J 

(10-4) 


no 


and  linearise,  disregarding  squares  and  higher  powers  of  the 
small  quantities  6w  and  c  .  The  result  is 

•  *  * 

A,qm1  +  (C-A)fl2m2  +  Q  c  1 3  -  Q2  c2  3  +  h^  -  h2  *  Li  , 

(10-6) 

*  %  * 

Anm2  -  (C  -  A)n2mi  ♦  ac23  +  a2 03 3  +  h2  +  nht  »  l2  . 


This  is  a  slight  and  strai ghtforward  generalization  of  eq.  (3-14) 
of  (Moritz,  1980b,  p.  17)  due  to  the  inclusion  of  the  relative 
angular  momentum,  which  is  necessary  since  it  is  not  convenient 
here  to  use  Tisserand  axes  for  the  entire  earth;  cf.  also  (Munk 
and  Macdonald,  1960,  p.  37). 

The  relative  angular  momentum  is  given  by 


jj  *  * 

earth 


hi 

h*> 


(10-7) 


This  is  eq.  (3-4)  of  (Moritz,  1980b,  p.  13),  but  with  ^  now 
replaced  by  ji  ,  which  is  necessary  since  the  relative  velocity 
was  denoted  by  u  in  the  earlier  report,  but  now  u  designates 
the  displacement  so  that  the  relative  velocity  is  u  ;  the  sym¬ 
bol  dM  denotes  the  mass  element. 

With  u  -  (u,v,w)  we  have 


1 

«< 

*  • 

1 

M 

<  • 

I 

r 

fi 

X  *  u  * 

• 

zu  -  xw 

. 

f  2 

xv  -  yu 
.  - 

_  M 

(10-8) 


Ill 


so  that 

hi  *  l\l  ftdM  •  (10’9) 

•arth 

We  split  up  this  integral  into  an  integral  over  the  core  and 
an  integral  over  the  mantle: 

ht  »  h l  -r  h®  (10-10) 

where 

hi  •  III  fidM  •  (10-11) 

core 

h®  »  |  j  f^M  .  (10-12) 

mantle 

Let  us  now  choose  the  coordinate  axis  to  be  Tisserand  axes 
for  the  mantle.  Tisserand  axes  are  defined  by  the  condition  that 
the  relative  angular  momentum  h  is  zero  (Moritz,  1980b,  p.  15). 
For  Tisserand  axes  for  the  mantle,  the  relative  angular  momentum 
due  to  the  mantle  only,  is  zero  : 

h®  *  0  .  (10-13) 


Thus 
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jj  ,iiM  * 

cor« 


(10-14) 


the  relative  angular  momentum  is  produced  by  the  core  only. 

The  substitution  of  (3-12)  into  (10-8)  gives  immediately 

fi  >  -<j  (v  -  e)xysinot  ♦  o  (v  -  ft  )y2  cosot  -  oftz2  cosot  (10-15) 

for  the  liquid  core.  We  substitute  into  (10-14)  and  perform  the 
integration  over  the  core.  We  integrate  over  the  unperturbed 
state,  for  which  the  core  is  an  ellipsoid  of  revolution,  the  z- 
axis  coinciding  with  the  axis  of  symmetry.  Then 

jj!  xydM  -  0  .  (10-16) 

core 

this  product  of  inertia  being  zero  because  the  z-axis  is  a  principal 
axis  of  inertia  (Heiskanen  and  Moritz,  1967,  p.  62).  We  further 
have 


C  * 
c 


* 

. 

core 


(x2  +  y2)dM  -  2 


» 

rf 

Jj 

y2dM 


core 


(10-17) 


because  of  rotational  symmetry;  Cc  is  the  polar  principal  moment 
of  inertia  of  the  core.  For  the  equatorial  principal  moment  of 
inertia  of  the  core  we  have 


A 

c 


III  (y2  +  z2)<im 

core 


C 

c 


f  r 


f 


core 


Z2dM  . 


(10-18) 
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Hence 


cor* 


(  10-19) 


Thus  the  substitution  of  (10-15)  into  (10-14)  and  Integra 
tion  gives,  in  view  of  (10-19), 


hi  ■  [£  o(v  -  8)  Cc  -  o$(Ac  -  |  Cc)]cosot 

or,  after  some  algebra  using  (2-58)  , 

hi«[fl(B-e)Ac  +  7  -  Ac)]cosot  , 

h2«[n(8-e)A  +  iov(C  -  A  )  ]  si  no  t  . 
c  c  c  c 


(10-20) 


(10-21) 


We  have  added  the  analogous  expression  for  h2  ,  which  can  be 
obtained  In  the  same  way. 

Now  we  have  to  determine  the  (negative)  products  of  inertia 
C13  and  c23  according  to  (10-5);  by  definition  (Moritz,  1980b, 
p.  9)  there  is 


(ft  f 

ci  3  •  -Dxz  ■  -  JJJ  xzdM  ,  C23  a  -Dyr  »  -  jjj  yzdM  .  (10-22) 

We  shall  calculate  them  separately  for  tidal  deformation  and  for 
rotational  deformation. 

Tidal  deformation.  Consider  the  v.ph«:  r  \  I  -ha  rittwn  i  <.  »1  *-  v «'  I  <»  (i 

ment  of  the  external  gravitational  potential  of  the  (deformed)  earth 
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•  n  ^.(cose) 

V(r.e,x)  •  V  l  —2— -  CA  eo*«x  ♦  B^slnmx)  .  (10-23) 

n-o  a*o  r 

By  (Heiskanen  and  Moritz,  1967,  p.  61,  «q.  (2-44c))  we  have 

A2  i  ■  GDx2  ■  “Gci3  »  B2 1  *  GDy*  ■  “SC2  3  ,  (  10-24) 

denoting  the  gravitational  constant  by  G  . 

The  gravitational  potential  of  the  undeformed  earth,  con¬ 
sidered  an  ellipsoid  of  revolution,  has  Dxs  «  0  »  0  .  Thus 

the  products  of  inertia  are  solely  due  to  the  deformation,  that 
is,  to  the  potential  Vi  ,  which  is  the  difference  of  the  gravi¬ 
tational  potentials  of  the  deformed  and  the  undeformed  earth.  For 
r  -  a  »  1  we  have  by  (8-18),  (8-19),  (9-8),  and  (9-9) 

Vx  «  kV#  «  kS  ■  k*s1necos9cos(o t  -  X )  (10-25) 


or 


Vi  ■  -j  kKP2i  (cose ) (cosotcos*  +  sinotsinx)  .  (10-26) 

We  put  r  ■  1  in  (10-23)  and  compare  the  terms  with  n  »  2  , 
m  «  1  with  (10-26).  This  gives  immediately 

Azi  *  y  ktecosot  ,  821  ■  y  kiesinot  ,  (10-27) 

so  that,  by  (10-24) 

c13  «-y  G"l»ekcosot  ,  c 2 3  »  -y  G'Uksinot  ,  (10-28) 

which  determines  cl3  and  C23  for  the  tidal  deformation. 

Rotational  deformation.  This  is  due  to  the  deformation  of 


the  earth  by  the  centrifugal  force,  similar  to  the  distorting 


115 


force  acting  on  a  rotating  unbalanced  wheel  (Moritz,  1980b,  p.  Zx). 
The  incremental  centrifugal  potential  by  (1*31)  is 

<t>  =  -Q2e  (xzcosct  +  yzsinot) 

=  -a2e  r2sinecos0cos(ct  -  a  ) 

On  the  earth's  surface,  for  r  =  a  ,  this  becomes 

$  =  -n2ea2s i ne cosecos (o t  -  a )  .  (10*29) 

The  comparison  with  (8-13)  and  (3-14)  shows  that  <t>  acts  like 
so.ie  kind  of  fictitious  external  potential  with 

k  =  -n2e  ;  (10-30) 

this  is  not  surprising  since  inertial  forces  such  as  the  centri¬ 
fugal  force  behave,  in  certain  respects,  like  genuine  forces.  Thus 
the  products  of  inertia  due  to  rotational  deformation  are  obtained 
by  substituting  (10-30)  in  (10-28)  : 

cR  ■  j  G”1  kn2ecosot  ,  C23  =  y  G" 1 kft2  e  s i  no t  .  (10-31) 

The  total  products  of  inertia  needed  in  (10-6)  are  the  sum 

T  D 

of  (10-28)  for  tidal  deformation  (they  should  be  written  C13 
and  cR)  and  (10-31)  for  rotational  deformation: 

C 1 3  *  cff  +  C?3  =  ~j  G_1k(K  -  n2e)C0Sat  , 

C23  *  ^2^  "**  cR  ~  -■j  G“ 1  k (k  -  fj2e)sinat 


(10-32) 


116 


Torque  components  L;  and  l_2  .  By  (10-25)  and  (10-26) 
we  write  the  tidal  potential  in  the  form 

Ve  a  3  *  P2  l  ( cose  )  (coso  tcosX  +  si  no  tsi  nX  )  .  (10-33) 

This  represents  V21  in  (Moritz,  1980b,  p.  3).  By  (ibid.,  p.  30) 
we  therefore  have 


L‘  ■  3MJ*  P21'{co-S9)  •e<i-90°>  • 

(10-34) 

■  ~3MJ*  P2T(ke7  v^°°)  • 

We  put  x  =  90°  and  x  =  0°  ,  respecti vely ,  in  (10-33)  and  sub 
stitute  in  (10-34)  .  The  result  is 


L 1  =  tcMJ2sinot  , 

L 2  =  -tcMJ2cosot  . 


(10-35) 


Euler-Liouvi  He  oondition.  Now  we  substitute  (10-21),  (10-32), 
and  (10-35),  together  with  (2-5),  into  (10-6).  We  get 


fl2Ce  -  (o  +  fl)[fiAe  +  nAc(0-e)  +  o(Cc~Ac)v  +  ■j£j3G"1ke  - 

(10-36) 


1 

7 


flG" liek  ]\  1 

si  not  ] 

\  =  kMJ2  1 

i 

s  i  not 

J 

-cosot  1 

l  ) 

^  -cosot 

By  (Heiskanen  and  Moritz,  1967,  p.  63)  we  have  (there  is  B  «  A 
because  of  rotational  symmetry): 
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,  ,  C  -  A 

02  =  J20  «  -pfpr  . 

whence 

kMJ2  *  <(C  -  A)  ( 10-37) 

since  a  *  1  .  Thus  (10-36)  gives 

n2Ce  -  (a  +  n)[nAe  +  QAC(6  -  e)  +  ^  a  (Cc  -  Ac)v  + 

+  6~lKk]  =  < ( C  -  A)  . 

After  some  straightforward  algebra,  we  finally  find 


e 


o  +  n 
Q 


fA  -  A 
_ c 

.  c 


8  + 


n2k 

v  +  -  e  - 

3GC 


K 

3GC 


C  -  A 

cn2 


(10-38) 


This  equation  relates  the  polar  motion  parameter  e  ,  the  core 
parameter  B  ,  and  the  Love  number  k  to  the  parameter  k  , 
characterizing  the  external  (tidal)  potential  by  (10-33);  note 
that  v  is  expressed  in  terms  of  B  and  e  by  (2-58)  : 


v 


2 


q  +  n 
0 


(10-39) 


It  is  a  consequence  of  Euler’s  and  Liouville's  equations  and  will, 
therefore,  be  called  Eul er-Li ouvi  11  e  condition. 
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It  is  sometimes  convenient  to  introduce  the  so-called 
"secular  Love  number"  kg  by 

k  «  3 G{ 1  0  96  (dimensionless)  (10-40) 

®  oi 

(Munk  and  Macdonald,  1960,  p.  26;  Moritz,  1980b,  p.  21).  Then, 
with  a  *  1  ,  eq.  (10-38)  takes  the  form 


a  +  a  A  -  A 

c 

£2  C 


+ 


C  -  A  k  ' 
C  k 

s 


C  -  A 
Cn2 


k  0  +  £2  1 


(10-41) 


For  tidal  frequencies  we  have  o  =  -£2  so  that  o  +  £2  will 
be  small;  so  is  ( C  -  A ) C  .  If  we  neglect  the  product  of  two  small 
quantities,  then  the  last  term  between  parentheses  on  the  left- 
hand  side  may  be  omitted;  whence 


o  +  £2 


e  - 


A-A  A  oC-A 

c  .  c  „  .  c  c 

- e  +  —  e  + - V 

C  C  2S2  C 


C  -  A 


0  +  S2 


£2 


(10-42) 


This  equation  is  identical  to  Molodensky's  (1961)  equation  (67). 


-rrit  r  I?  rrfrrl,‘*ijavii'  - 


a»  •tariff 


Equation  (10-33),  or  its  approximate  form  (10-42),  repre¬ 
sents  the  seventh  condition  needed  to  uniquely  determine  the  7 
parameters  h ,k  ,  1  ,<*  ,8  ,y  ,e  as  mentioned  at  the  end  of 
sec.  9. 

Determination  of  parameters.  The  relevant  equations  are 
thus  (9-33),  (9-34),  and  (10-3S).  We  can  solve  them  as  follows. 

We  first  take  from  the  system  (9-33)  the  equations 

M(b )  =  0  , 

p-i  N(b)  +  R(b)  +  H (b ) W ' (b )  *  b2o (o +2n )< “ 1  $  =  0  ,  (10-43) 

core 

2bR(b )  -  P(b)  +  3bW' (b)H(b)  -  0  . 

Since  by  (8-25)  the  quantities  H(b),  M(b),  N(b),  P(b),  and 
R ( b )  are  linear  combinations  of  the  Love  numbers  h,  k,  1  with 
constant  coefficients,  the  system  (10-43)  has  the  form 

aih  +  a2k  +  33!  *  a0  , 

b  1  h  +  b2k  +  b3l  *  bo  +  b4o(o+2n)B  ,  (10-44) 

Cih  +  c2k  +  c3l  =  c0  • 

We  can  solve  this  equation  for  h,  k,  1,  giving  them  as  linear 

functions  of  o(a+2n)8  . 

Mol odensky  (1961)  has  considered  two  models,  both  more  real¬ 
istic  than  our  simple  earth  model  featuring  a  homogeneous  liquid 
core  and  an  elastic  mantle.  Molodensky's  two  models  have  a  more 
detailed  core  structure:  Model  1  has  a  liquid  core  whose  density 
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increases  towards  the  center,  and  Model  2  has,  in  addition,  a 
solid  (elastic)  inner  core.  Still,  Molodensky's  Model  1  is 
sufficiently  similar  to  our  simple  model,  so  that  his  numerical 
values  can  be  used  as  an  illustration  of  the  present  formulas. 
Molodensky  (1961)  gives  for  Model  1  : 

h  *  0.6206  +  0.4711  *  10‘3  ?.l?  0  , 

k  *  0.3070  +  0.2384  *  10“3  0  3  ,  (10-45) 

1  =  0.0904  -  0.0112  x  10-3  0  . 


These  formulas  show  clearly  that,  for  a  liquid  core  model,  the 
Love  numbers  depend  on  frequency  (the  strong  dependence  of  0 
the  frequency  o  is  shown  in  the  table  below). 

Now  we  can  use  the  second  equation  of  (9-33)  together  with 
(9-34),  (10-38),  and  (10-39).  We  have 


e2b*”lv  -  H(b)  *  0  , 
\  (y  -  0  )e2  + 


*  '■  n '  ->2  4.  £Jt_ £  g  -  £  e  =  0  , 


+  n 

A  -  C 

c 

fc  -  A  0  C  ) 

C  C  J.  c 

n 

C 

.  C  0  Cj 

C  -  A 

K  ^ 

m 

<  ■  -  -  ; 

3GC 

_ 

c  n2 

n2k 

0  +  -  c 

3GC 


(10-46  ) 


on 


the  last  equation  is  obtained  by  eliminating  v  from  (10-38) 
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by  means  of  (10-39).  Since  H(b)  is  a  linear  combination  of 
h,  k,  1  and  since  h,  k,  1  are  functions  of  s  by  (10-45), 
(10-46)  represents  a  system  of  three  equation'  for  8  ,y  ,e  . 

If  it  has  been  solved,  we  finally  get  a  from  the  third  equation 
of  (9-23). 

For  Model  1  with  a  surface  flattening  of  1/298.3  and  the 
core  values 

e2  *  0.005092  ,  Cc/C  *  0.  i  180  .  (  10-47  ) 

Molodensky  gives  the  following  expressions: 

41.87 

6  =  -  +  1.9  ,  (10-48) 

0.2136  -  100(o  +  R)/o 


C  q  a  +  R 

—  «  1  +  —2—  (8  -  4.3)  -  .  (10-49) 

e  0  C  -  A  R 

Eq .  (10-49)  is  to  be  understood  in  the  following  way.  Denote 
by  e o  the  value  of  e  for  a  rigid  earth  (solid  throughout). 
Then  k  *  0  because  of  rigidity,  and  Ac  55  0  *  Cc  because  the 
core  is  missing.  Thus,  for  a  rigid  earth,  (10-38)  reduces  to 


o  +  R  A 

IT"  rj  e° 


C  -  A 
=  *Tr2' 


(10-50  ) 


This  equation  determines  c0  .  The  quantity  q  in  (10-49)  has 
the  definition 
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.  A2 a  a  centrifugal  force  at  equator  f  10-51  i 

q  g  gravity  at  equator  *  '  ' 

it  reduces  to  a*  for  units  in  which  a  ■  1  ,  g  •  1  . 

Thus  e/e o  represents  the  ratio  between  polar  motion  ampli¬ 
tudes  for  a  liquid  core-elastic  mantle  model  and  for  a  rigid  earth. 

Numerical  values  can  be  found  in  the  following  table,  taken 
from  Molodensky  (1961)  for  Model  1. 


o  +  O 

a 

6 

e/e  o 

Precession 

0 

198 

1 

Principal  nutation 

-1/6800 

212 

0.996 

(18.7  years ) 

+1/6800 

185 

1.003 

Semiannual  nutation 

-1/183 

-125 

00 

CO 

o 

•-H 

+1/183 

57 

! 

1.036 

Fortnightly  nutation 

-1/13.7 

-4 

1.080 

+1/13.7 

+  7 

1.026 

This  table  can  be  compared  with  the  table  on  p.  122  of  (Moritz, 


1980  b) 


123 


for  the  Poincarfe  model  (rigid  shell  and  liquid  core);  our  present 
e /e o  corresponds  to  U/Uo  there.  The  present  values  are  much 
more  in  agreement  with  observation.  The  comparison  shows  that  ela¬ 
sticity  of  the  shell  reduces  the  effect  of  the  liquid  core  as 
compared  to  a  rigid  shell;  this  fact  has  already  been  pointed  out 
by  Jeffreys  (  1949 ) . 

Note  the  strong  dependence  of  the  (dimensionless)  parameter 
B  on  frequency.  Therefore,  B  is  sometimes  called  Molodensky)s 
resonance  parameter  for  the  core. 

The  proper  frequencies  correspond  to  free  oscillation  without 
external  force;  for  them  we  have  <  •  0  .  Indeed,  from  V@  *  0 
there  follows  k  ■  0  by  (8-16).  As  in  the  Poincarfc  model,  there 
are  four  proper  frequencies:  the  axial  spin  mode  (ASM),  o0  ,  the 
Chandler  wobble  (CW),  o!  ,  the  nearly  diurnal  free  wobble  ( NDFW ) , 
a 2  ,  and  the  turn-over  mode  (TOM),  a3  ;  cf.  (Smith,  1977;  Moritz, 

1980b,  p.  116).  The  "trivial"  modes  ASM  and  TOM  are  independent  of 
the  constitution  of  the  body;  they  are  also  now  the  same  as  for  a 
rigid  body,  namely 


o0»0,  o  3  ■  - n  ;  ( 10-52 ) 

cf.  (Moritz,  1980b,  p.  83). 

The  proper  frequency  for  NDFW  is  obtained  by  putting  *c  =  U 
in  (10-42).  Molodensky  gives  for  his  Model  1  the  value 

a 2  a  -1.00214O  (Molodensky  1).  (10-53) 


This  may  be  compared  to  the  value 


o2 


-  n 


1  + 


T  THT 


(10-54) 


ib&iZ* 111* 


-  asaiM 
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for  the  Poincarfc  model,  cf.  (Moritz,  1980b,  p.  116),  with 

e  *  e2/2  and  C  «  C-C  .  With  the  values  (10-47)  this  gives 

n  c 

o 2  -  -1.00289  fl  (Poincari)  .  (10-55) 

For  the  Chandler  frequency  (CW)  we  cannot  neglect  o+a 
so  that  we  must  use  (10-41).  In  fact,  the  Chandler  frequency  is 
so  small  that  we  can  neglect  it  as  compared  to  n  ,  so  that  now 

•  1  ♦  •  1  .  (10-56) 

The  second  equation  of  (10-46)  may  be  written 

j  (>  -  »)  •*§■  ♦  »  -  «  ■  0  .  (10-57) 

For  small  o/n  ,  the  first  term  may  be  neglected,  and  there 
remains 

»  ;  «T7  •  '  (l0'58> 

From  (2-58)  there  follows 

Jjj  -  6-e*e-e»0  .  (10-59) 

With  <  »  0  (proper  frequency,  no  external  potential)  and  (10-58) 
and  (10-59),  eq.  (10-41)  reduces  to 


iH Mi 
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5+a  f  A-A{  Aq  n  C  -  A  k 
n  l  C  C  £j  +  o C  C  ka 


We  divide  by  c  end  solve  for  o  .  The  result  is,  writing  o  *oj 
for  CW  : 


n 


(10-60) 


To  an  error  of  about  1  X  ,  this  equals 

oi  C-A  w 

«—  -  -  1  -  JL. 

n  A  -  A  k 

c  s 


(10-61) 


a  result  given  by  Molodensky  (1961).  For  his  Model  1  he  gets  a 

value  of  Oi  which  corresponds  to  a  Chandler  period  of  433  days. 

Introducing  the  Euler  frequency  o_  by 

£ 

°e  *  n  (10-62) 


(Moritz,  1980b,  p.  10),  (10-60)  may  be  written 


(10-63) 


126 


Relation  to  McClure's  theory,  This  theory  (McClure,  1973), 
also  described  In  (Moritz,  1980b),  Is  based  on  an  entirely  elastic 
model  without  liquid  core.  Such  a  model  can  be  considered  a  limit 
of  a  Molodensky  model  for  vanishing  core: 

b  -  0  ,  Ac  -  Cc  ■  0  .  (10-64) 

Then  (10-63)  for  the  Chandler  frequency  reduces  to 


1  - 


oi 


k 

*7 


1  ♦ 


He 

1C 


. - -  a 

is.  E 

a 


(10-65) 


identical  to  eq. 
Eq.  (10-41) 


(4-5)  of  (Moritz,  1980b). 
reduces  to 


e 


o  n 
Q 


-*■ 


C  -  A 


k 

TT 

s 


C  -  A 
*  "IF" 


l 


ifiij  .(10-66) 


We  multiply  by 


C 

* 


(10-67) 


(by  (10-62)) and  introduce  the  tidal  frequency  and  the  nuta- 

tional  frequency  Au^  by 


UK  ■  -0  , 

Aw i  ■  a  -  n  *  - ( o  +  fl ) 


(10-68) 
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(Moritz,  1980b,  pp.  95  and  97).  We  get 


1 

c 

<J,  Aw.  1 

l  ♦  JE.  +  _JL 

f  k  °E 

1  + - i 

«- 

1 

°E 

a  <  - 

r  k  Aw.  ) 

n  + - i 

a  0 

m 

l  a 

a3 

c; 

ft 

(10-69) 


Using  (10-65),  this  Is  readily  seen  to  be  the  same  as 


c 


E 


+  °C) 


(10-70) 


This  equation  Is  equivalent  to  the  tidal  term  In  eq.  (5-24) 
of  (Moritz,  1980b,  p.  34),  as  we  shall  see  by  relating  the  con¬ 
stant  <c  to  the  expansion  of  the  tidal  potential  (ibid.,  sec.  1). 
For  the  tesseral  part  V21  of  degree  2  and  order  1  we  have 
(ibid.,  eqs  .( 1-10)  and  (5-14)): 


v2l  «  2a2P21  (cose)  J  B^sin^v  +  8  ^  +  x) 


(10-71) 


the  sum  being  extended  over  all  tidal  frequencies  w^  ,  and 
and  denoting  amplitude  and  phase,  respectively. 

On  the  other  hand,  in  the  present  report  we  have  used  the 
form  (2-43)  for  the  tidal  potential: 


*  «r2sinecos9cos(ot -  x) 

«  -JK  r2P2I(cose)cos(ot  -  x) 


(10-72) 
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The  comparison  with  (10-71)  shows  that  (10-72)  represents  any 
term  of  the  sum  (10-71)  with  o  ■  -u  by  (10-68)  and 

k  ,  S3  «  90°  (10-73) 

with  r  «  a  »  1  on  the  earth's  surface. 

Now,  McClure's  theory,  eq.  (5-24)  of  (Moritz,  1980b),  limited 
to  the  tidal ly  induced  term  with  frequency  ,  gives 


By  (1-28)  we  have 

m  ■  mj  +  im2  *  c(cosat +  i sinot)  ■  eel0t  .  (10-75) 

By  (10-73)  and  (10-75)  we  thus  see  that  (10-70)  is  indeed  equi¬ 
valent  to  (10-74),  as  was  to  be  shown. 

Relation  to  Poincarfe's  theory.  Let  us  specialize  eqs.  (10-46) 
for  the  Poincarfe  model  featuring  a  rigid  mantle  and  a  homogeneous 
liquid  core.  Since  the  rigid  mantle  is  not  deformed,  we  have 

H(b)  •  0  ,  (10-76) 


and  the  first  equation  of  (10-46)  gives 
y  *  0 


(10-77) 
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The  second  equation  then  yields 

SJlSL  g  -  —  e  «  i  e2$  ,  (10-78) 

0  o  c 

or 

8  +  |  (6  -  e)  ■  |  e2g  ,  (10-79) 

identical  to  (4-64 ).' Since  the  rigid  mantle  is  not  deformed,  the 
gravitational  potential  does  not  change,  so  that 

V:  »  0  (10-80) 

and  therefore 
k  *  0  . 

Hence  the  third  equation  of  (10-46;  becomes 

a  +  n  '  A  -  C  fC-Ao  C  1  I  C-A 

- - - —  e  +  — - -  -  +  —  B  =  < -  . 

a  C  C  n  C  Cn2 

w  — 

In  sec.  4  we  have  shown  that  (10-79),  or  (4-64),  is  equivalent 
to 

Acoe  +  (Aco  +  Ccfl)(8  -  e )  =  0  ;  (10-83) 

cf,  (4-67).  Using  (4-61)  and(4-69),  namely 

u0  *  ,  v0  *  0(8  -  e )  ,  ( 10-84 ) 

this  has  been  brought  to  the  form  (4-54), 

Ac<ju0  +  ( Aca  +  Ccn)v0  *  0  ,  (10-85) 

which  is  one  of  the  two  basic  equations  of  Poincar&'s  theory. 

Let  us  now  show  that  the  second  basic  equation,  namely  (4-53), 
can  be  derived  from  the  Euler-Li ouvi 1 1 e  condition  (10-82).  We  put 
«  *  0  for  free  oscillations.  Then,  by  elementary  manipulations. 


(10-81) 


(10-82  ) 
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we  can  bring  (10-82)  into  the  form 

[Ao  -  (C  -  A)0  -  (o  +0)Cc)  +  !  Ac  +  (C  -  AJ  s  =  0  . 

o  +  Q  L  &  -1 

(10-86  ) 

Since  a  +  n  is  very  small  for  tidal  frequencies,  we  may  write 

(a  +  0)Cc  *  (a  +  fl)Ac 
and  neglect 

Sr1  (C=-A=> 

as  the  product  of  two  small  quantities.  Thus  (10-86)  becomes 

[Ao  -  (C-A)O  -  A  (a  +  ft )  ]  e  +  A  (o  +  fl)B  •«  0  .  (10-87) 

C  v 

This  can  be  written  in  the  form 

[Ao  -  (C-A)a]e  +  A o(o  +  n)(8-e)  =  0  .  (10-88) 

With  (10-84)  this  finally  becomes 

[  Ao  -  (C-A)n]u0  +  Ac(o  +  n)vo  =  0  ,  (10-89) 

identical  to  (4-53).  This  completes  the  derivation  of  the  fundamen¬ 
tal  equations  of  Poincar&'s  theory  from  Molodensky's  theory,  showing 
the  equivalence  of  both  (completely  different!)  theories  for  the 
Poi ncar§  model . 


,-!;  ."*  .t,V.'  -  ‘  PI 
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11 .  Recent  Developments 


Molodensky's  theory  has  been  serving  as  a  starting  point  for 
a  number  of  subsequent  Investigations.  It  has  been  refined  and  made 
more  general,  but  also  simplified  and  made  more  elegant.  Thus  Molo- 
densky  was  pioneering  in  this  fie>d,  as  he  was  in  the  geodetic 
boundary-value  problem. 

The  work  of  Shen  and  Manslnha  (1976).  They  provided  a  uni¬ 
form  treatment  of  mantle  and  core  by  using  the  equations  of  elasti¬ 
city  (cf.  sec.  6)  throughout  the  earth,  regarding  the  liquid  core 
(or  liquid  outer  core)  as  a  limiting  case  for  an  elastic  core  with 

v  *  0  (11-1) 

(see  end  of  sec.  5).  This  requires  the  consideration  of  the  ellip- 
ticity  in  the  elastic  equations,  since  the  ellipticity  is  essential 
for  core  resonance  as  we  have  repeatedly  seen.  The  motion  of  the 
liquid  core  is  treated  in  terms  of  spheroidal  and  toroidal  oscilla¬ 
tions  in  the  way  outlined  in  sec.  3. 

Shen  and  Mansinha's  results  agree  well  with  Molodensky's. 
Naturally,  these  authors  have  been  able  to  use  recent  better  models 
for  the  earth's  interior,  especially  the  core;  cf.  also  (Melchior, 
1980a; Smi th  ,  1980). 

The  method  of  Sasao,  Okubo  and  Saito  (1980).  They  gave  a 
particularly  elegant  formulation  of  Molodensky's  by  generalizing  the 
well-known  Poincarfe  equations  for  a  rigid  mantle  and  a  liquid  core; 
cf.  (Moritz,  1980b),  eqs.  (12-22)  with  (12-32): 
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•  • 

Au  -  i(C-A)ftu  +  Ac(v  +  iflv)  *  L  , 

(11-2) 

A  u  +  A  v  +  iC  Qv  *  0 

c  c  c 

Here  A  and  C  are  the  principal  moments  of  inertia  of  the  whole 
earth,  Ac  and  Cc  those  of  the  core,  the  complex  number  L 
combines  the  and  x2  components  of  the  lunisolar  torque: 


L  =  L  i  +  i  1.2  ,  (11-3) 

and  the  complex  numbers  u  and  v  represent  the  rotation  of  the 
.  earth  and  that  of  the  core,  respectively  (ibid.  ,  p.  118);  see  also 
sec.  4  of  the  present  report. 

(Sasao  et  al . ,  1980)  succeeded  in  finding  a  surprisingly 
simple  equivalent  of  eqs.  (11-2)  for  Molodensky's  problem.  Using 
the  same  notations  as  above,  it  reads  (ibid. ,  eqs.  (37)  and  (32)): 

•  •  • 

Au  -  i(C-A)flu  +  Ac(v+iflv)  +  n(c+inc)  ■  L  , 

(11-4) 

A  u  +  A  v  +  iC  nv  +  nc  *  0 
c  c  c  c 

The  new  feature  is  the  appearance  of  the  complex  numbers 
c  «  Ci 3  +  ic23  » 

(11-5) 

cc  -  cf 3  +  ic?3  , 

combining  the  products  of  inertia  C;3  and  c2a  for  the  earth 


.  ■  ... 
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(cf.  Moritz,  1980b,  pp.  16-17)  and  for  the  core,  respectively. 

These  quantities  can  be  expressed  in  terms  of  the  functions 
(8-25)  characterizing  the  elastic  behavior  of  the  mantle,  by 
means  of  the  boundary  conditions  discussed  in  sec.  9. 

In  fact,  the  first  equation  of  (11-4)  is  equivalent  to  the 
Euler-Liouvi 1 le  condition,  (10-38)  or  the  third  equation  of  (10-46). 
whereas  the  second  equation  of  (11-14)  is  a  more  elegant  version 
of  the  second  equation  of  (10-46);  cf.  the  relation  to  Poincar&'s 
theory  discussed  in  sec.  10. 

The  theory  of  (Sasao  et  al.,  1980),  as  well  as  that  of  (Shen 
and  Mansinha,  1976),  Is  applicable  to  more  general  earth  models 
than  those  used  by  Molodensky;  again,  the  numerical  results  are 
similar. 

The  work  of  Smith  and  Wahr.  Smith  (1977)  used  the  normal 
modes  (proper  modes,  proper  frequencies)  of  the  earth  to  study 
polar  motion  and  nutation.  The  most  extensive,  detailed  and  accurate 
application  of  this  theory  was  made  by  Wahr  (1979,  1981a, b,c). 

Wahr  (1981a)  used  differential  equations  of  motion  for  an 
elastic  earth  (6-33)  in  the  form  (6-37).  For  the  liquid  outer 
core  the  same  equations  hold  with  (11-1),  so  that  a  uniform  treat¬ 
ment  of  the  whole  earth  is  possible.  Also,  the  ellipticity  is  con¬ 
sistently  taken  into  account. 

The  differential  equations  of  motion  have  the  form 

•*  • 

p  (u^  +  x  u)  3  _L jj  ,  (11-6) 

where  1L  is  a  linear  operator  acting  on  the  displacement  vector 
jj  involving  first  and  second  derivatives  with  respect  to  ^ 

The  coordinate  system  used  is  the  uniformly  rotating  nutation 
frame;  cf.  secs.  1  and  6. 
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These  equations  are  solved  for  subject  to  the  boundary 
conditions  (cf.  secs.  8  and  9)  that  the  normal  displacement  u  •  n 
(n  is  the  unit  vector  normal  to  the  boundary  surface),  the  stress 
vector  (5-13),  the  change  of  gravitational  potential,  V i  , 

and  the  quantity  (9-24)  must  be  continuous  across  any  boundary. 
Furthermore,  the  displacement  vector  as  such  must  be  continu¬ 
ous  across  any  "welded  boundary".  Such  a  welded  boundary  separates 
elastic  regions  of  different  density,  etc.,  such  as  crust  and 
mantle;  the  core-mantle  boundary  is  not  a  welded  boundary  since 
the  core  consists  of  a  liquid  whose  movement  is  restricted  only 
by  continuity  of  the  normal  displacement  (it  can  more  freely  along 
the  boundary ) . 

Eqs.  (11-6)  are  transformed  to  spherical  coordinates,  with 
ellipsoidal  corrections  being  taken  into  account.  The  solution  u 
is  consequently  expressed  in  the  form  (3-1)  : 


H  -  II(i“  ♦  £)  .  (11-7) 

no 

as  a  linear  combination  of  spheroidal  and  toroidal  modes  (again, 
with  small  corrections  for  ell ipticity) . 

For  nutation,  the  terms  with  m  *  1  are  responsible.  Wahr 
uses  the  truncated  form 

u  -  l\  +  s\  +  T3  (11-8) 

for  his  purpose.  This  form  is  also  employed  by  Shen  and  Mansinha 
(1976),  whereas  Molodensky's  theory  is  equivalent  to  using 

u  «  li  +  S2  ;  (11-9) 


22^22 
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The  toroidal  term  t}  may  be  interpreted  as  a  small  spatial 
rotation.  In  fact,  (3-10)  can  be  written  as 

ij  «  9  X  X  (11-10) 

with 

Tr" 1 s i not 

Tr^cosot  ’»  (11-11) 

0 

this  representation  holds  generally  for  any  t|  .  Such  a  vector 
product  is  characteri sti c  for  an  infinitesimal  rotation;  cf. 
(Moritz,  1980b,  p.  78).  If  the  vector  £  does  not  depend  on  the 
radius  vector  r  ,  then  we  have  a  rigid  rotation  of  the  body  as 
a  whole.  This  is  the  case  if 

Tr”1  »  const.  or  T(r)  *  cr  ,  (11-12) 

c  being  a  constant.  (Otherwise,  the  spheres  r  *  const,  under¬ 
go  different  rotations.) 

For  a  homogeneous  liquid  core,  T  is  indeed  proportional 
to  r  ,  as  (3-17)  shows.  A  rigid  mantle  obviously  also  rotates 
as  a  whole,  so  that  (11-12)  holds.  Of  course,  the  rotations  of 
core  and  mantle  will  in  general  be  different,  so  that  the  constant 
c  differs  for  core  and  mantle. 
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Remarkably  enough,  the  relation  (11-12)  also  holds,  at  least 
approximately,  for  far  more  general  earth  models  featuring,  e.g., 
an  elastic  mantle,  a  liquid  outer  core,  and  an  elastic  inner  core. 
Each  of  these  three  regions  is  undergoing  a  nearly  rigid  rotation; 
cf.  Fig.  7  in  (Smith,  1977). 

Using  (11-10)  we  may  write  (11-8)  in  the  form 

u  *  0  x  x  +  u'  (11-13) 

The  first  term,  t|  ,  expresses  a  small  rotation  which  will,  by 
definition,  be  called  nutation  ;  the  remainder 

u'  «  +  l\  (11-14) 

will,  again  by  definition,  be  called  the  effect  of  body  tide. 

This  splitting-up  of  the  total  displacement  into  nutation  and  tide 
is,  to  a  certain  degree,  arbitrary  but  is  in  agreement  with  pre¬ 
vious  definitions  of  polar  motion  and  nutation  as  we  shall  see 
i  iow.  A  similar  decomposition  has  been  considered  in  sec.  3;  cf. 
(3-13)  and  (3-21). 

Thus  nutation  (rotation)  and  tidal  deformation  have  been 
separated.  Due  to  nutation,  the  earth's  surface  (again  assumed 
spherical)  rotates  like  a  rigid  sphere.  Hence,  as  far  as  geometry 
(kinematics)  is  concerned,  the  rigid  body  theory  of  Smith  (1977) 
and  Wahr  (1979,  sec.  7),  as  extensively  discussed  in  (Moritz,  1980b, 
'  i .  .  nd  11),  may  again  be  applied.  In  particular,  the  formulas 
for  precession  and  polar  motion  ( ibid ,  ,  eqs.  (11-37)  and  (11-39)) 
can  be  used  as  far  as  they  refer  to  pu  ely  kinematically  defined 
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axes.  This  is  true  for  the  rotation  axis,  so  that  the  first 
equation  of  (11-37),  ibid.,  remains  valid: 

nR  »  ijj-  9  ,  (11-15) 

where 

9  *  9 i  +  i 02  , 
n  *  r\i  +  in2  , 

K 

1 

2  *  *R  "  ®3  * 

|  n3 

Thus  the  vector  q  represents  the  very  small  difference  between 

R 

e„  ,  the  unit  vector  of  the  rotation  xis,  and  e3  ,  the  coordi- 
nate  unit  vector  of  the  z-axis  of  the  reference  frame  (the  nutation 
frame  as  mentioned  above)  which  has  a  constant  direction  in  space. 
The  complex  number  nR  combines  the  first  two  components  of  the 
spatial  vector  n^  (the  third  component,  n3  ,  is  easily  seen  to 
be  negl igibly  smal 1 ) . 

With  the  figure  axis  F  we  must  be  careful.  From  (Moritz, 
1980b)  we  take  the  second  equation  of  (11-37): 

np  .-10  ,  (11-19) 


(11-16) 

(11-17) 

(11-18) 


which  is  equivalent  to 
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Sr  ‘  *  *3  '  1  x  *3  i  (11-20) 

cf.  ibid.,  eqs,  (11-15)  and  (11-35).  Thus  ep  is  obtained  by 
rotating  the  nutation  axis  03  by  6_  ;  cf.  (11-10)  with  £3 
taken  for  x.  •  In  other  terms,  e^,  is  the  unit  vector  of  the 
z-axis  to  which  the  body  tide  ju '  refers.  Assuming  for  the  mo¬ 
ment  that  u/  refers  to  the  body  frame  as  defined  at  the  end 
of  sec.  6  and  used  by  Molodensky,  we  can  identify  the  axis  F 
with  the  z-axis  of  the  body  frame;  more  about  this  will  be  said 
below.  In  particular,  the  axis  F  is  an  average  axis  of  figure 
rather  than  the  instantaneous  axis  of  maximum  inertia. 

Polar  motion  is  referred  to  the  body  axis  e^.  ,  hence  it  is 
expressed  by  the  complex  number 

p*  ■ ",  -  "r  ■ (  H1*  •  (n-21) 

in  agreement  with  ibid. .  eq.  (11-39). 

The  analogy  with  rigid-body  motion  breaks  down  for  the  angular 
momentum  axis  H  ,  which  is  not  a  kinematical  but  a  dynamical 
quantity.  Hence,  (11-37),  ibid,  ,  does  not  hold  for  nH  .  However, 
since  the  basic  relation  between  the  angular  momentum  vector  H 
and  the  torque  L  , 

$  ■  t  •  <n-22> 

is  independent  of  the  physical  constitution  of  the  body,  the  rela¬ 
tion  ( 13-37 ) , i bi d .  ,  valid  for  a  rigid  body  as  well  as  for  Poincar&'s 
model  , 


must  hold  also  In  the  present  case.  The  correspondi ng  polar  motion 
is  then  given  by 


with  (11-19)  and  (11-23).  j 

This  behavior  of  the  kinematical  axes  R  and  F  and  of  the 
dynamical  axis  H  is  fully  analogous  to  the  case  of  the  Poincarfe  j 

model;  cf.  (Moritz,  1980b,  pp.  128  -  130). 

Relation  to  Kinoshita's  theory.  Throughout  the  present  re 
port  we  have  restricted  our  considerations  to  one  individual 

frequency 


only;  the  general  case  is  then  a  sum  of  the  individual  frequency 
components ;  cf.  also  (Moritz.  1980b.  pp.  95  -  98  and  140  -  141). 
Thus,  e.g.,  e  has  the  usual  exponential  form 


where  the  amplitude  n  is  a  real  constant  if  (11-12)  holds  (or 
a  function  of  r  otherwise);  for  generality  we  have  added  a 

phase  angle  y  . 

This  complex  quantity  8  holds  for  the  Mol odensky-Wahr  model. 
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For  the  rigid  earth  we  have  a  similar  quantity  given  by  (ibid. t 
eq .  (11-33))  and  denoted  by  er  ;  its  amplitude  will  be  desig¬ 
nated  by  nr  .  Wahr  (1979,  1981c)  writes 


and  uses  his  theory  only  to  calculate  the  small  quantity  (n-nr)/nr; 
nr  itself  is  obtained  from  the  extremely  accurate  rigid-body 
theory  of  Kinoshita  (1977);  cf.  also  (Moritz,  1980b,  sec.  9). 

The  computational  advantages  of  this  procedure  with  regard  to 
accuracy  are  evident. 

Remark  on  reference  systems.  The  body  frames  used  in  (Moritz, 
1980b)  and  in  the  present  report  have  always  been  Tisserand  axes. 

For  rigid  body,  such  axes  are  fixed  to  that  body,  so  that  a  rigid 
body  is  at  rest  with  respect  to  Tisserand  axes.  For  a  nonrigid 
body,  Tisserand  axes  are  defined  in  such  a  way  that  the  body  is 
"on  the  average"  at  rest  with  respect  to  such  a  system. 

A  plausible  mathematical  formulation  of  "being  at  rest  on 
the  average"  is  obtained  by  condition 


(H-27) 


minimizing  the  residual  distortions  u'  in  (11-13);  p  is  the 
density  and  dv  is  the  volume  element.  This  condition  can  be 
shown  to  be  basically  identical  to  (ibid. ,  eq .  (3-7)). 

The  obvious  choice  of  the  volume  of  integration  v  is 
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v  ■  the  whole  earth  .  (ll-28a) 

This  choice  Is  appropriate  for  a  completely  elastic  earth  with¬ 
out  liquid  core;  the  body  axes  xyz  In  secs.  3  to  8  of  (Morit2, 
1980b)  have  been  defined  In  such  a  way. 

In  the  case  of  a  liquid  core  it  is  better  to  link  the  axes 
to  the  mantle  only,  by  taking 

v  ■  the  mantle  (plus  crust)  .  (ll-28b) 

The  symmetry  axes  of  the  rigid  mantle  used  as  coordinate  frame 
for  Poincar&'s  theory  ( ibid. ,  secs.  12  -  13)  are  such  axes,  and 
so  is  the  body  frame  used  in  Molodensky's  theory:  the  condition 
(10-13)  is  equivalent  to  (11-27)  with  ( 1 1  - 2  8  b ) ;  cf.  (ibid. .  p.  15). 
A  third  choice  is 

v  «  the  crust  (ll-28c) 

(which  may  practically  be  identified  with  the  earth's  surface). 

In  view  of  plate  motion,  etcs.,  it  may  be  argued  which  of  the  de¬ 
finitions  (ll-28b)  or  ( 1 1 -28c )  is  geophysically  more  appropriate; 
for  an  excellent  discussion  see  (Smith,  1981). 

At  any  rate,  for  a  spherically  symmetric  elastic  mantle,  both 
definitions  (11-2 8b )  and  (ll-28c)  are  equivalent  provided  the 
condition  (11-12)  holds  since  then  the  mantle  (plus  crust)  rotates 
as  a  whole. 

For  Wahr's  theory,  the  definition  (ll-28c)  is  most  natural 
since  we  have  seen  that  it  amounts  to  interpreting  tJ  a 
rigid  rotation  of  the  spherical  earth's  surface.  If  according  to 


— -^'T  -i  r 
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(11-13)  we  split  off  the  effect  of  this  rotation,  £  *  x  , 
from  the  displacement  u_  ,  the  residual  displacement  u/  will 
refer  to  a  system  which  contains  no  such  rotation,  that  is,  no 
term  t|  ,  and  this  may  be  shown  also  to  hold  for  Tisserand 
axes  (another  expression  of  the  body's  "being  on  the  average  at 
rest"  in  such  a  system).  Thus  u'  refers  to  "crustal  Tisserand 
axes"  (ll-28c),  which  practically  coincide  with  Molodensky's  t 
"mantle  Tisserand  axes"  since  (11-12)  holds  to  very  high  accuracy- 
Hence  Wahr's  u_'  is  practically  identical  to  Molodensky's  ui  , 
both  referring  to  a  body  frame  which  is  a  crust-mantle  Tisserand 
frame,  whereas  Wahr's  original  u.  refers  to  the  nutation  frame 
as  we  have  pointed  out  above. 

In  the  notation  of  eq.  (11-20),  the  unit  vector  eF  has  the 
direction  of  the  2-axis  of  the  body  frame  thus  defined.  As  a 
Tisserand  axis  it  is  optimally  stable  with  respect  to  the  body, 
whereas  the  real  instantaneous  figure  axis,  the  instantaneous  axis 
of  maximum  inertia,  is  very  unstable.  This  situation  is  completely 
analogous  to  the  purely  elastic  case  (Moritz,  1980b,  pp.  40  and 
49  -  50).  As  a  matter  of  fact,  the  z-axis,  represented  by  the 
vector  ,  can  be  considered  an  average  axis  of  figure.  This 
is  particularly  obvious  from  Fig.  6.1.  ( ibid . ,  p.  38):  the  daily 
average  of  the  instantaneous  figure  axis  F  is  Fo  ,  and  the 
average  of  F0  over  a  Chandler  period  is  0  ,  corresponding  to 
the  z-axis. 

For  the  sake  of  completeness  we  point  out  that  also  the 
instantaneous  rotation  axis  R  represents  an  average,  although 
in  a  completely  different  sense.  In  fact,  in  a  nonrigid  body, 
each  volume  element  undergoes  a  differential  rotation  described 
by  the  tensor  R  given  by  (3-23)  or  (5-5),  which  in  general  varies 
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from  point  to  point.  Such  a  differential  rotation  is  superimposed 
on  the  general  body  rotation  and  makes  the  rotational  speed  slight 
ly  different  at  each  point. 

Thus  even  the  "instantaneous  rotation  vector  w  of  the  earth 
as  a  whole"  can  only  be  an  average  over  the  whole  earth  of  the 
rotation  vectors  of  all  volume  elements.  Jeffreys  (1970,  sec. 

7.04)  proposes  to  define  w  by  the  condition 


minimum 


(11-29) 


where  is  the  velocity  of  any  point  with  respect  to  an  inertial 
system.  This  definition  is  intimately  related  to  the  definition 
of  Tisserand  axes,  and  our  pR  (11-21)  corresponds  precisely  to 
it. 

Thus  the  rotation  axis  R  represents  a  space  average  but  is 
instantaneous  in  time,  whereas  the  body  axis  z  mentioned  above 
is  average  of  the  figure  axis  over  a  long  time  period. 
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PART  C 


JEFFREYS'  VARIATIONAL  APPROACH 


12 .  Variational  Principles 


The  variational  approach,  initiated  by  Jeffreys  (1949), 
has  led  to  the  poineering  work  (Jeffreys  and  Vicente,  1957a, b), 
which  gave  the  first  numerical  results  for  nutation  and  tidal 
Love  numbers  for  an  earth  model  consisting  of  a  liquid  core  and 
an  elastic  mantle.  To  my  knowledge,  this  method  has  never  been 
used  again,  partly  because  these  papers  are  very  difficult  to 
read,  as  has  been  remarked  by  several  authors,  beginning  with 
Molodensky  (1961).  Besides  its  historical  importance,  however, 
Jeffreys'  method  has  considerable  theoretical  interest,  and  to 
the  patient  and  detective-minded  reader,  the  brilliant  obscurity 
of  his  papers  is  a  source  of  fascination  (though  sometimes  mixed 
with  frustration) . 

We  shall  here  try  to  provide  a  simple  introduction  to  Jeffreys 
approach,  pointing  out  the  basic  ideas  and,  perhaps,  encouraging 
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and  facilitating  the  study  of  the  original  papers.  We  begin  with 
some  well-known  facts  from  theoretical  mechanics,  partly  follow¬ 
ing  (Courant  and  Hilbert,  1953,  pp.  242-249). 

Hamilton's  principle.  Consider  a  mechanical  system  of  n 
degrees  of  freedom,  described  by  n  parameters  (generalized  or 
lagrangian  coordinates)  qlf  q2,  ...»  qn  .  Then  the  motion  of 
the  system  between  two  instants  t0  and  ti  of  time  t  is 
such  that 


J  (T  -  U)dt  ■  extremum  .  (12-1) 

to 

Here  U  denotes  the  potential  energy,  considered  as  a  function 
of  the  coordinates  and  the  time* 


U  «  U(q1#  q2,  ....  qn,  t)  *  U(qr,  t)  ,  (12-2) 

and  T  denotes  the  kinetic  energy,  which  is  quadratic  in  the 
velocities  qr  *  dqr/dt  : 


a  q  q 
rs^r^s 


(12-3) 


where,  in  general,  the  coefficients  a  are  functions  of  the 

JT  8 

coordinates  and  the  time  : 


(12-4) 
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In  (12-3)  and  in  the  following  we  shall  use  the  summation 
convention,  summing  over  repeated  indices  (here,  r  and  s  ) 
from  1  to  n  . 

Introducing  the  Lagrangian  function  L  by 
L  »  T  -  U  »  L(qr,  qr,  t)  ,  (12-5) 

we  may  write  (12-1)  as 


t  * 

Ldt  ■  extremum  (12-6) 

to 


or 


5L  -  0  ,  (12-7) 

an  extremum  being  equivalent  to  a  vanishing  variation  6L  as 
known  from  variational  calculus.  This  condition  leads  to  Lagrange 1 s 
equations  : 


d  aL  at 
it  a<5^  ~  aqr 


(12-8) 


An  important  special  case  is  that  of  small  variations.  In  this 
case,  the  a  in  (12-3)  are  constants,  and  the  potential  energy 
(12-2)  is  a  quadratic  function  of  qr  .  Thus 
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T  *  7  ar.V.  V 
U  *  7  br.V.  * 


(12-9) 


with  constant  coefficients  ars  and  bfs  .  With  ( 12-9 )t  Lagrange ' s 
equations  (12-8)  give 

*„V  ■>„<(,  ■  0  >  (12-10) 

a  system  of  n  linear  second-order  ordinary  differential  equations 
with  constant  coefficients  for  the  unknown  coordinate  functions 
qr  *  qr(t)  . 

Vibrating  string.  This  case  provides  a  continuous  analogue 
to  the  discrete  case  of  n  degrees  of  freedom  (Fig.  12.1).  The 


u 


string  is  kept  fixed  at  the  two  end  points  A  and  B  .  Then  the 
kinetic  energy  T  and  the  potential  energy  U  are  given  by  the 
expressions 


F!83B?5! 


i 


? » 


$-  i 
£  i 


t 

p 


where 


u  * 


i 

(12-12)  1 


p  and  u  are  constants  denoting  density  and  tension  of  the 
spring,  respectively,  and  i  •  AB  . 

The  expressions  (12-11)  may  be  considered  continuous  analogues 
of  (12-9),  the  integrals  corresponding  to  the  sums  in  (12-9) 
(implied  by  the  summation  convention). 

Hamilton's  principle 


Ldt  * 
co 


(T  -  U)dt  *  extremum 

< 

to 


(12-13) 


with  (12-11),  or  <5L  *  0  ,  can  be  shown  by  standard  methods  of 
variational  calculus  to  lead  to  the  linear  partial  differential 
equati on 


p U  -  yu  * 
tt  xx 


(12-14) 


i 

* 


v 


describing  the  motion  of  a  vibrating  string,  which  may  be  solved 
for  the  displacement 
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u  «  u(x,t)  .  (12-15) 

Vibrating  membrane.  The  step  from  the  one-dimensional  string 
to  the  two-dimensional  membrane  is  evident:  we  have 


T 

U 


1 

7 


1 

7 


'  f  * 

pu2dxdy  , 

i 

G 

u(u2  +  u2)dxdy  , 
G 


(12-16) 


G  being  the  area  of  the  membrane.  Hamilton's  principle  <SL  * 

«  6  (T  -  U )  =*  0  leads  to  the  linear  partial  differential  equation 

»utt  -  “<“**♦  V  - 0  <12-17> 

for  the  displacement 

u  ■  u(x.y.t)  .  (12-18) 


Jeffreys'  variational  principle.  Now  we  go  to  three-dimensions 
and  a  displacement  vector  ^  ■  (u1.U2.U3)  .  Then  the  kinetic  energy 
is  expressed  by 


T 


1 

7 


f  p | u | 2dv 


1 

7 


p  (u2  +  u2  + u|)dv 


(12-19) 
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where 

dv  *  dxdydz  ■  dx^x^dxj  (12-20) 

is  the  volume  element  and  v  the  volume  of  integration  (in  our 
case,  the  earth).  The  potential  energy  is  given  by  the  expression 


92  V  9  ( V  +V  )  9  V  9  U  a 

-  U  .  U.  +  - - - —  ut - U  .  - -  + 

9Xi9X  1  2  9X±  9Xi  1  9  X 


9  V 


9  X 


9  U 


9X  , 


9U  . 

- 1  _  ^ 

P  - 

3  9Xj 


dv 


(12-21) 


where  V  denotes  the  gravitational  potential  and  p  the  density 
as  used  in  sec.  6;  the  stress  tensor  pi;j  is  expressed  in  terms 
of  u±  in  the  usual  way  (6-36). 

The  expression  (12-19)  for  the  kinetic  energy  is  immediately 
obvious,  being  a  natural  generalization  of  the  first  equation 
of  (12-16).  The  quantity  0  is  quadratic  in  £  and  its  derivatives, 
similarly  to  the  second  equation  of  (12-16),  but  is  otherwise 
much  more  complicated.  Jeffreys  and  Vicente  (1957a)  derived  it 
from  energy  considerations.  We  shall  here  be  satisfied  with  the 
fact  that  it  provides  the  correct  differential  equations. 

In  fact,  the  condition  6L  *  0  leads  by  standard  variational 
methods  to 
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3U  3V  3u.  3  V  3 2  V 

PU  «  -p  -  -  +  p - +  0  u.  —  • —  + 

3Xk  3Xi  3XA  3Xk  3Xi3Xk 

3  3Pik 

+  P  -  (V  +  SI.)  +  — ,  (12-22) 

3xi  3xk 

which  is  the  system  of  equations  of  motion  for  an  elastic  non¬ 
rotating  earth,  identical  to  (6-33)  together  with  (6-37),  for 
u  »  0  and  W  «  V  ,  correspondi ng  to  the  absence  of  rotation  and 
centrifugal  fo. ce.  This  system  is  to  be  solved  for  the  displace¬ 
ment  vector 

u i  *  Ui(x1  ,x2  ,X3  ,t)  .  (12-23) 

In  this  way  the  earth  is  considered  a  three-dimensional  vibrating 
body  whose  oscillations  are  excited  by  the  lunisolar  force  (the 
free  oscillations  of  the  earth  are  also  described  by  these  equations, 
with  »  0  ) .  The  analogy  to  the  vibrating  membrane  is  obvious. 

Finally  we  point  out  an  important  feature  of  the  variational 
method  :  the  condition  <5L  -  0  gives  not  only  a  differential 
equation  such  as  (12-22),  but  also  boundary  conditions  (continu¬ 
ity  of  the  normal  displacement,  of  the  stress  vector,  etc.).  These 
boundary  conditions  are  here  obtained  in  a  formal  mathematical 
way,  without  recourse  to  physical  considerations  such  as  used  in 
secs.  8  and  9  . 
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13.  Method  of  Jeffreys  and  Vicente 


Jeffreys  (1949)  and  (Jeffreys  and  Vicente,  1957a,  b)  used 
the  variational  approach  in  an  original  and  interesting  way.  They 
did  not  use  the  principle  6L  ■  0  to  derive  (12-22)  from  (12-19) 
and  (12-21),  although  they  mention  this  possibility,  but  they  try 
to  reduce  the  continuous  to  the  discrete  case  at  the  very  begin¬ 
ning.  We  shall  here  try  to  outline  the  basic  principles,  referring 
the  reader  to  the  original  paper  for  further  details. 

The  idea.  Starting  again  from  the  principle 


Ldt 

to 


r:l 

(T  -  U)dt  ■  extremum 
to 


(13-1) 


we  try  to  reduce  L  from  the  continuous  case,  for  which  (12-19) 
and  (12-21)  hold,  to  the  discrete  case  (12-9).  This  is  done  in 
the  following  way. 

Assume  that  the  elastic  displacement  vector  u.  can  be  re¬ 
presented  as  a  linear  function  of  a  finite  number  n  of  variables 

(Lagrangian  parameters)  qi,q2 . qn.  Th,>s  is  possible,  e.g., 

by  using  the  representation  (11-7)  in  a  truncated  form  such  as 
(11-8).  We  thus  write 


u . 


y  c,  q  *  c.  q 
L,  ir^r  ir^r 


1 


(13-2) 


using  the  summation  convention:  indices  such  as  i  or  j  will 
run  from  1  to  3  ,  corresponding  to  the  coordinates  in  three- 
dimensional  space,  whereas  r  or  s  run  from  1  to  n  .  Let 
us  make  explicit  the  functional  dependence: 

ut  ■  u1(x1 ,x2  * x3 » t )  , 

Cir  *  Cir(Xi,X2,x3)  ,  (13-3) 

qr  «  qr(t)  . 

The  coefficients  cir  are  considered  known,  the  parameters  qr 
unknown . 

Let  us  illustrate  these  general  relations  by  the  concrete 
case  (11-8).  The  A-component  of  t{  ,  by  (3-5)  with  (11-12),  has 
the  form 

(Tih  *  *  T(r)cosasin(ot  -  A) 

*  -rcosecosA  •  c  sinot 

+  rcosesinA  •  c cosot  .  (13-4) 

The  first  term  on  the  right-hand  side  then  contributes  to  the 
sum  (13-2)  with 

c.  ■  c.  (xi,x2,x3)  -  rcosesinx  , 

JL  i  a.  j. 

(13-5) 

*  qr(t)  *  “ c sinot 
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Thus  qr  represents  essentially  the  unknown  coefficient  c  .  (It 
is  clear  that  spherical  coordinates  r,e,x  can  be  expressed  in 
terms  of  xx.x2.x3,  so  that  cir  is  indeed  a  function  of  the 
latter;  also,  there  is  hardly  any  danger  of  confusing  the  sub¬ 
script  r  with  the  radius  vector.)  Similar  contributions  come 
from  the  second  term  of  (13-4)  and,  in  fact,  from  any  component 
of  S°  and  Tm  in  (11-8). 

Having  thus  clarified  the  nature  of  cir  and  qr  ,  we  substi¬ 
tute  (13-2)  into  (12-19)  and  (12-21).  The  differentiations  3/3xi, 
etc.,  are  performed  on  the  known  functions  cir  ,  and  then  the 
integrations  over  the  earth's  volume  y  are  carried  out,  which 
do  not  affect  the  functions  qr(t)  because  they  are  independent 
of  position  xi  .  The  result  is  a  Lagrangian  L  *  T-U  of  form 
(12-5)  which  is  quadratic  in  qr  because  the  original  Lagrangian 
was  quadratic  in  ui  .  Then  Lagrangian's  equations  (12-8)  give 
a  system  of  ordinary  linear  differential  equations  for  qr  , 
similar  to  (12-10) . 

Consideration  of  rotation.  This  barest  sketch  of  the  basic 
idea  must  now  be  made  more  precise  and  more  concrete.  First  of  all, 
we  must  consider  a  rotating  frame  of  reference.  In  fact,  the  initial 
equations  (12-19)  and  (12-21)  refer  to  a  nonrotating  inertial  frame. 
They  are  transformed  to  the  uniformly  rotating  nutation  frame  which 
Jeffreys  and  Vicente  are  using.  Furthermore  we  must  admit  small 
rotations  of  the  earth  with  respect  to  this  frame,  described  by 
the  vector  £  as  in  (11-13)  or,  equivalently,  by  small  angles 
I'.n'.m'  as  in  (4-8).  We  thus  get  an  expression  of  form 


L  «  T  -  U 


I  •  •  •  » 


1  '  ,m,'  n '  )dv  , 


(13-6) 


1S5 


similar  to  the  one  obtained  by  subtracting  (12-21)  from  (12-19) 
but,  In  addition,  containing  rotation.  The  displacements  ui 
refer  to  the  body  frame  and  are  thus  the  same  as  ut  in  secs. 

7  to  10,  as  well  as  the  same  as  u^  In  (11-13);  Jeffreys  and 
Vicente  write  uj^  for  our  ui  in  (13-6). 

Secondly,  we  split  up  the  Integration  over  the  earth  Into 
Integrations  over  the  core  and  over  the  mantle: 


v  cor*  mantle 


The  mantle.  Eqs.  (7-17)  give  for  the  earth's  surface 
r  ■  a  «  1  : 

ur(l)  »  hS  «  hV#  ,  (13-8) 

ue(l)  «  13S/36  «  l3Ve/39  ,  (13-9) 

here  we  have  used  (8-10),  (8-11),  and  (8-16),  putting 

ve  "  Va>  *  ved>  •  (13-10) 

Here  V  is  taken  in  the  form  (10-72)  which  depends  on  the 

© 

tidal  frequency 


o  -  -w  .  (13-11) 

The  Love  numbers  also  depend  on  frequency,  so  that  (13-8)  should 


I 


i 


_ 
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(13-17) 
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where 


qi(t)  =  -j  n2  I  hj  Bj  sin(a)3 1  +  b3  )  , 
q2(t)  =  -j  n2 1  hjBjcos(o)j  t+ ej )  . 

In  the  same  way  we  treat  (13-9),  obtaining 


(13-18) 


ue  ( 1 )  -  q3(t)  ♦  q4(t)  (13-19) 

wi  th 

q  3  ( t )  1  jBjsin(«jt  +  e3 )  , 

j  (13-20) 

q4(t)  *  -j  n2  l  1  jBjcos(uijt  +  B3 )  . 

Consider  finally  the  gravitational  potential  Vi  induced 
by  the  tidal  potential  Ve  .  By  (8-18),  (8-19),  and  (13-10)  we 
have 


V i ( 1 )  »  kVe  .  (13-21) 

Again,  this  is  valid  for  one  frequency  only,  and  the  total  poten¬ 
tial  Vi  is  the  sum  of  all  such  frequency  contributions.  In  the 
usual  way  we  get 
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V  i  ( 1 )  *  ftn(e»A)q5(t)  +  Su(8,x)q6(t) 

with 

q5(t)  *  ~  Q2  l  k  jBjsin(Wjt  +  , 

qs(t)  *  |  n2  I  k 3B  .cos(Ujt  +  e3)  • 


(13-22) 


(13-23) 


The  functions  qi(t),  ...,q6(t)  will  be  taken  as  Lagrangian 
parameters  for  the  mantle.  For  the  reader  familiar  with  spectral 
analysis  it  is  obvious  that  they  are  essentially  the  transforms 
of  the  Love  numbers 


hj  *  h(Uj)  ,  kJ  *  k(Uj)  ,  ^  «  1 (wj)  (13-24) 

from  the  frequency  domain  into  the  time  domain. 

Next  we  consider  the  displacement  vector  u  for  general 
r  s  1  .  On  expressing  the  Love  numbers  by  (13-8),  (13-9),  and 
(13-21)  for  one  particular  frequency  : 

h  *  u  r  ( 1 )  /  S  , 

1  *  u0(l)/S0  .  (13-25) 

k  =  Vx (1)/S  , 

substituting  into  (8-25),  and  inserting  the  resulting  functions 
H(r)  and  L(r)  into  (7-17),  we  find  ur,ug,u^  as  linear 


functions  of  ur(l),u0(l),  and  V i  ( 1 )  .  If  we  perform  this 
operation  for  all  frequencies  and  sum  the  results,  we  get 

in  view  of  (13-17),  (13-19),  and  (13-22)  expressions  of  the 
form 

ur(r,9,X,t)  =  Q?(r,e ,X )qg(t)  ,  (13-26) 

summation  over  s  from  1  to  6  being  implied,  and  similar¬ 
ly  for  u0  and  ux 

Finally  we  obtain  the  rectangular  components  U1.U2.U3  by 
the  rotation  (3-8),  which  gives 

ui(xi,x2,x3,t)  *  Uis (xi ,x2 ,x3 )qg (t)  ,  (13-27) 

expressing  the  displacements  as  linear  functions  of  the  Lagrangian 
variables  qs  ;  the  functions  Uia  are  known. 

Our  choice  of  qg  ,  based  on  ur(l),u0(l),  and  Ve(l)  , 
has  been  motivated  by  the  wish  to  have  a  close  relation  to  the 
Love  numbers  h,k,l  so  that  the  developments  in  secs.  7  and  8 
could  be  used.  Jeffreys  and  Vicente  (1957a)  use  a  slightly  differ¬ 
ent  choice,  based  on  ur ( 1 )  ,ur  (b) ,  and  Ve(l)  ,  ur(b)  being 
the  radial  displacement  at  the  core-mantle  boundary. 

The  core.  By  (4-4)  we  have 

ut  =  aijXj  (13-28) 

(we  may  write  x^  instead  of  x°  since  u^  is  small),  where 
ai;j  are  given  by  (4-13)  in  the  nutation  frame.  As  usual  we  put 
n  *  0  *  n',  obtaining  Uj.  as  linear  functions  of  l,m,T,m'. 
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In  the  body  frame  which  we  are  using  now,  the  even 

depend  on  1  and  m  only  (the  last  rotation  (4-8)  is  then 
missing).  Thus 

ui  “  Ai]  +  Bim  *  '13-29) 

the  coefficients  being  known  functions  of  position  as  usual 
(even  linear  in  xi )  . 

Lagrange's  equations.  We  now  substitute  (13-27)  for  the 
mantle  and  (13-29)  for  the  core  into  (13-6)  and  perform  the 
integration.  The  result  has  the  form 

L  =  L(di  »d2i  ■  •  •  »<Uo  )  (13-30) 


with  the  Lagrangian  variables  qi  to  q6  as  above  and  with 


qv(t)  -  1* 
q9(t)  -  m' 
q9(t)  «  1 


(13-31) 


q  i  o  ( t )  =  m 


L  will  be  quadratic  in  qr  for  the  reasons  outlined  at  the 
beginning  of  this  section;  more  precisely  it  is  found  to  have 
the  form 

1  .  .  1 

L  »  —  a  oo-ib  qo  +  c  oa  -  da 
2  rs^r^s  2  rs^r^s  rs^r^s  rMr 


(13-32) 
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with  constant  coefficients  ars*-**»dr  •  how,  of  course,  the 
summation  indices  r  and  s  each  run  from  1  to  10  . 

The  lagrangian  equations  (12-8)  give  immediately 

•t,’,  +  <e.t  ■  +  *  dr  ’  0  •  (13'33> 

The  further  treatment  is  standard.  It  is  best  to  use  complex 
combinations  (note  sin  and  cos  in  expressions  such  as 
(13-18) l)1 } 

qi  +  iQa  a  M 

q3  +  iq4  =  X2 

Qs  +  ic)6  55  X3 

q7  +  iq8  *  x4 

q9  +  iqio  *  X5 


> 


> 

(13-34) 

* 


% 


« 


and  seek  the  solution  in  the  form 

xk  *  x°eiot  (k  =  1,2 . 5)  (13-35) 

(in  view  of  (13-11)  this  amounts  to  a  transformation  into  the 
frequency  domain;  a  phase  factor  is  without  importance). 

Then 


1)  The  deeper  reason  for  the  possibility  of  complex  combination 
Is  the  rotational  symmetry  of  the  earth  model  used. 
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ar  «  ioXk  ,  Xk  «  -o2Xk  ,  (13-36) 

and  (13-33)  reduces  to  a  system  of  linear  algebraic  equations. 

By  (13-18)  we  have  for  one  frequency  o  »  -w  only  : 

Xi  «  io2hjBje"i(“;Jt+ (13-37) 

so  that,  apart  fom  constants,  x°  essentially  is  h  (we  again 
write  h  instead  of  h^  or  h(w^)  ).  Similarly,  x£  and  x° 
essentially  are  the  Love  numbers  1  and  k  ,  respectively.  Our 
present  x4  and  x5  are  immediately  seen  to  be  identical  to  x' 
and  x  in  (4-29)  and  (4-46),  so  that  (4-42)  and  (4-58)  show 
that  x°  and  x°  are  essentially  (apart  from  costant  factors) 
identical  to  Molodensky's  parameters  e  and  6 

The  system  of  linear  algebraic  equations  to  which  (13-33) 
reduces,  thus  has  the  form 

A i j h  +  A 1 2 ^  +  A i 3 1  +  Ai4e  +  A i 5 B  *  Bj  , 

A2  1  h  +  A2 2 k  +  A23I  +  A2 4 e  +  A 2 5 S  *  B2  , 

(13-38) 

A  5 1 h  +  A  5  2  k  +  A  5  3 1  +  A  5  4  e  +  A55S  *  B  5  , 

It  is  comparable  to  the  system  of  5  equations  obtained  by  elimi 
nating  y  from  the  system  of  6  equations  (10-43)  and  (10-46) 
for  h,k,l,e,B,y  and  should  be  equivalent  to  it  as  a  whole  (al- 


though  the  individual  equations  are  different). 

Each  of  the  equations  of  sec.  10  has  been  derived  separately 
from  physical  considerations  (boundary  conditions  and  Euler- 
Llouville  equation).  The  present  system  ( 13- 38),  however,  has  been 
obtained  as  a  whole  directly  from  the  variational  principle 
6 L  *  0  using  (13-32),  in  a  formal  way  without  physical  considera¬ 
tions.  (This  is  reminiscent  of  the  fact  that,  also  in  the  contin¬ 
uous  case,  the  condition  <$L  «  0  gives  not  only  the  respective 
differential  equation  but  also  correspondi ng  boundary  conditions, 
see  the  remark  at  the  end  of  sec.  121) 

This  difference  between  Moi uuensky ' s  method  and  the  method 
of  Jeffreys  and  Vicente  implies  both  advantages  and  drawbacks. 

In  Molodensky's  method,  each  condition  has  a  simple  physical  mean¬ 
ing,  it  is  thus  more  transparent  and  can  be  modified  and  extended 
to  more  complex  earth  models  (e.g.,  containing  a  liquid  outer  core 
and  an  elastic  inner  core).  Jeffreys'  approach  has  a  conceptually 
simpler  structure;  once  the  Lagrangian  (13-6)  has  been  established, 
everything  else  follows  in  a  logically  strai ghtforward  manner. 
Nevertheless,  the  details  are  enormously  complicated,  and  only 
an  author  of  the  physical  insight  and  mathematical  skill  of  Sir 
Harold  Jeffreys  could  devise  such  an  approach  and  lead  it  to  a 
successful  conclusion. 

In  view  of  these  difficulties,  Jeffreys  and  Vicente  (1957a, b) 
consider  only  two  simplified  earth  models:  the  central  particle 
model,  consisting  of  a  homogeneous  core  and  a  central  mass  point 
that  is  to  represent  the  solid  inner  core,  and  the  Roche  model 
using  a  continuous  density  distribution  in  the  core  according  to 
Roche's  law.  For  the  same  reason,  beginning  at  an  early  stage, 
computations  are  performed  numerically  instead  of  analytically: 
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for  instance,  already  the  coefficients  a.  . in  the 

Lagrangian  (13-32)  are  calculated  numerically.  Equations  such 
as  (13-38)  have  no  physical  interpretation.  Thus  this  approach, 
though  logically  more  elegant,  is  physically  less  transparent. 

Also,  as  Molodensky  (1961)  has  remarked,  it  is  difficult  to  get 
a  clear  idea  of  the  approximations  involved  in  the  present  method. 

The  approach  of  Jeffreys  and  Vicente  has  not  been  followed 
by  later  investigators  and  is  out  of  fashion  nowadays.  Still, 
the  variational  approach  to  partial  differential  equations  in 
general  is  quite  fashionable  in  both  pure  and  applied  mathematics 
because  it  is  theoretically  powerful  *nd  computationally  convenient. 
Thus  it  does  not  seem  impossible  that  Jeffreys'  approach  will  some 
day  be  revived  in  a  different  form. 


14 .  Numerical  Aspects 


Let  us  finally  give  some  numerical  values  showing  the 
practical  implication  of  the  theories  discussed  in  the  present 
report.  We  shall  use. the  following  abbreviations: 

JV1  (Jeffreys  and  Vicente,  1957a) 

JV2  (Jeffreys  and  Vicente,  1957b) 

Ml  (Molodensky,  1961),  Model  1 

M2  (Molodensky,  1961),  Model  2 

SM  (Shen  and  Mansinha,  1976) 

W  (Wahr ,  1979,  1981  b,c) 

K  (Kinoshita,  1977),  rigid  earth  . 

Nutati on .  Astronomers  use  the  angles  Ae  (nutation  in 
obliquity)  and  Ai|i  (nutation  in  longitude)  :  e  4  23.5°  is  the 
obliquity  of  the  ecliptic.  They  are  related  to  nR  or  np  by 
eq.  (11-55)  of  (Moritz,  1980b),  with  e  =  e  . 

From  (McCarthy  et  al . ,  1980),  Table  3,  we  take  the  follow 
ina  values  for  the  nutation  of  the  rotation  axis  n_  : 


18.6  year  Ac 
Arsine 


9 

6 


.2050"t0.0017 

,8409"±0.0025 


semiannual  Ac 
Arsine 


0.578"  ±0.004" 
0.533"  ±0.004" 


fortnightly  Ac 
Arsine 


0.0912"±0.0015 
0.0859 "±0.0010 


"  9.2273 

"  6.8713 


9.1972  9.2145 

6.8327  6.8556 


9.1951  9.1985 

6.8312  6.8356 


0.5500 

0.5046 


0.5706 

0.5212 


0.5377 

0.4859 


0.5710 

0.5216 


0.5686 

0.5196 


4 


0.0886 

0.0813 


0.0914 

0.0835 


0.0912 

0.0835 


0.0913 

0.0834 


0.0908 

0.0831 


Note  that  these  rigid-earth  values,  serving  as  a  reference, 
do  not  correspond  to  Kinoshita. 

(Wahr,  1981c,  p.  725)  gives  for  the  nutation  of  the  mean 
figure  axis,  np  (which  is  the  z-axis  of  the  body  system,  a 
Tisserand  axis  for  mantle  or  crust)  : 


Term 

18.6  year  Ac 
Arsine 


semiannual  Ac 
Ai^si  ne 


M2 


SM 


9.2278" 

6,8743" 


9.2044 

6.8441 


9.2012 

6.8400 


0.5534" 

0.5082" 


0.5719 

0.5232 


0.5745 

0.5253 


W 


9.2025 

6.8416 


0.5736 

0.5245 


0.0 

0.0 


0.0972 

0.0899 


fortnightly  Ae 
Arsine 


1949" 

1881" 


0.0978 

0.0904 


0.0977 

0.0905 
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We  see  that  It  Is  necessary  to  take  liquid-core  effects  into 
account  if  we  want  to  get  agreement  with  observations.  A  rigid- 
earth  theory  Is  no  longer  adequate.  On  the  other  hand,  various 
modern  liquid-core  theories  give  very  similar  results;  still, 
the  differences  on  the  order  of  0.002"  are  becoming  observation- 
ally  significant. 

At  Its  General  Assembly  in  Montreal  in  August  1979,  the 
International  Astronomical  Union  (IAU)  has  made  the  transition 
from  a  rigid  earth  to  a  liquid-core  model  and  adopted  nut.ational 
constants  based  on  the  Molodensky  2  model.  The  International 
Union  of  Geodesy  and  Geophysics  (IUGG),  at  its  General  Assembly 
in  Canberra  in  December  1979,  has  asked  the  IAU  to  reconsider 
this  decision.  It  now  appears  probable  that  the  IAU  will  reverse 
its  decision  and  adopt  values  based  on  a  Wahr  model. 

Nearly  diurnal  free  wobble.  Let  this  frequency  be  denoted 
by  c2  (cf.  (10-53));  then  the  following  values  for  o2/n  nave 
been  obtained: 


-o2/n 

JV1 

1.00224 

JV2 

1.00403 

Ml 

1.00214 

M2 

1.00216 

W 

1.00218 

Chandler  period.  Empirical  determinations  give  a  period 
of  about 
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( 435 ±3 )  sidereal  days 

(Dahlen,  1980).  Correspondi ng  values  from  earth  models  are  (in 
sidereal  days) 


di rect 

corrected 

for  ocean 

JV1 

4 

392 

430 

0  V  2 

395 

433 

Ml 

401 

433 

M2 

402 

436 

U 

403 

The  effect  of  the  oceans  on  the  Chandler  period  is  considerable. 
Dahlen  (1976)  found  that  oceans  act  to  increase  the  Chandler 
period  by  28  sidereal  days.  See  also  (Smith,  1977). 

Love  numbers.  They  are  frequency  dependent.  We  shall  give 
them  for  some  main  diurnal  tidal  frequencies.  For  the  notation 
Oj,  Pj,  Ki  cf.  (Melchior,  1978,  pp.  27  and  50-51);  they  cor¬ 
respond  to  nutaLional  periods  of  13.7  days  (fortnightly), 

183  days  (semiannual),  and  »  (precession),  respectively. 

Love  number  h 


J  V 1 

JV2 

Ml 

M2 

W 

0i 

0.584 

0.603  ' 

0.617 

0.614 

0.603 

Pi 

0.555 

0.568 

0.594 

0.593 

0.581 

<i 

0.492 

0.551 

0.527 

0.535 

0.520 

.T'itfi  ■  -Ti^  W-fia.  iatentAsjli 


xUIuIi 


1.081 
J.  082 
93  0.084 


These  values  have  been  taken  from  (Molodensky,  1961)  for  JV1, 

JV2 ,  Ml  and  M2,  and  from  (Wahr,  1979);  cf.  also  (Melchior,  1978, 
pp.  155  and  158). 

Concluding  remarks.  These  numerical  values  are  shown  only 
for  the  purpose  of  illustrating  numerical  differences  between 
different  theories.  Since  the  various  theories  should  be  essential¬ 
ly  equivalent,  the  discrepancies  will  be  due  primarily  to  different 
models  for  the  earth's  interior  and  to  computational  inaccuracies. 

The  earth  models  used  range  from  the  highly  simplified  models 
employed  by  Jeffreys  and  Vicente  to  modern  sophisticated  earth 
models  used  by  Wahr.  Therefore,  and  because  of  the  precision  of 
theory  and  computations,  Wahr's  results  should  be  most  reliable. 

A  discussion  of  the  physical  aspects  of  earth  models  used, 
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important  as  the  subject  is,  is  beyond  the  scope  of  the  present 
report,  and  the  reader  is  referred  to  the  literature,  e.g.  (Smith, 
1977;  Melchior,  1980a, b;  Wahr,  1981b). 

In  all  these  models,  the  earth  is  considered  a  sphere  or  an 
ellipsoid  of  revolution,  and  the  mass  distribution  is  to  have  a 
corresponding  symmetry.  Deviations  of  the  actual  earth  from  such 
an  idealization  can  cause  considerable  discrepancies.  For  example, 
we  have  already  mentioned  the  influence  of  the  oceans  on  the 
Chandler  period;  ocean  loading  and  other  local  effects  also  have 
a  considerable  influence  on  earth  tides.  The  basic  references 
for  a  geophysical  discussion  of  the  earth's  rotation  are  (Munk 
and  Macdonald,  1960)  and  (Lambeck,  1980). 
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